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Abstract 

This is the first part of lecture notes on my approach to field quantization. I explain on a simple 
scalar-field model the physical motivation and show some preliminary applications (field produced 
P-f by a pointlike charge and the vacuum-to- vacuum loop diagram). The approach is based on an 

appropriate construction of a representation of harmonic-oscillator Lie algebra. 
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Sensible mathematics involves neglecting a quantity when it is small — not neglecting it just 
because it is infinitely great and you do not want it! 
P. A. M. Dirac (1975) fl 

For many the infinities that still haunt physics cry for further and deeper quantization, but 
there has been little agreement on exactly what and how far to quantize. 
D. R. Finkelstein (2006) Q 



I. INTRODUCTION 



When solving problems in quantum mechanics one often notices that discrete energy levels 
are determined by practically one mathematical condition — finite norm of eigenvectors. The 
reason for the finiteness is simple — squared norm represents sums of probabilities associated 
with maximal sets of jointly measurable physical quantities. It should be normalizable to 1 
and thus cannot be infinite. Finiteness of another important quantity, (ip\H 2 \i(j), guarantees 
that fluctuations AH = \J (H 2 ) — (H) 2 are finite. From a mathematical point of view 
(i/j\H 2 \ifj) < oo means that and H\ip) belong to the same Hilbert space, a condition 
that defines the domain of H. As we can see, the apparently trivial "law of finiteness of 
quantities that should be finite" is raised in quantum mechanics to the level of important 
technical principle. It determines values of physical quantities (spectra of operators) and 
sets of acceptable states of physical systems (domains of operators). If one allowed for 
discrete spectra and infinite norms then harmonic oscillators would possess infinitely many 
arbitrarily large negative energy levels. 

In quantum field theory we become much more tolerant. Infrared and ultraviolet "catas- 
trophes" (infinities occurring at various stages of calculations) are regularized and renor- 
malized. Algorithms for extracting physically reasonable numbers from infinite theoreti- 
cal predictions were invented and resulted in several Nobel Prizes in physics. The proce- 
dures generally seem to work, but how to explain to a beginner why we cannot proceed 
in a similar way in standard quantum mechanics with "divergent probabilities associated 
with non-normalizable eigenvectors"? We could take the inverted Gaussian associated with 
E = —Hcu/2, regularize the divergent norm by a cutoff, divide the eigenvector by such a 
regularized norm, compute average energy, and finally remove the cutoff. The procedure 
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would certainly imply some kind of "catastrophe" , but basically what would be wrong with 
it? A similar level of logical rigor seems acceptable in quantum field theory. 

Yet, in spite of all the successes of quantum field theory I'm convinced that it is the 
paradigm of quantum mechanics that is physically correct. We have to determine mathe- 
matical structures by the well-definiteness of physical quantities. Not all mathematically well 
defined models are physical, but models that are mathematically inconsistent are unphysi- 
cal for sure. The fact that it is practically impossible to perform a nontrivial quantum field 
theoretic calculation without encountering an ill defined mathematical expression shows, in 
my opinion, that the most fundamental theory of physics is formulated in a wrong way. 

What I write is, of course, not original. Dirac repeatedly stressed that a deep change is 
needed in conceptual foundations of field quantization. In his opinion quantum field theory 
is still at the stage analogous to atomic physics before the advent of Schrodinger's equation 
- we had agreement between theory and experiment (Bohr's model of atom predicted 
correct spectra), but theoretical ideas were completely wrong. In his last two papers 3, J], 
published in the year of his death, he formulated his desiderata for fundamental physical 
theories. He believed that quantum dynamics should always be described by Heisenberg 
equation of motion with appropriate Hamiltonian. Modifications are expected in algebraic 
properties of dynamical variables since here freedom is immense. He stressed the unexplored 
potential of reducible representations of fundamental symmetries. 

In this context it should be reminded that the Heisenberg equation is completely unrelated 
to classical Euler-Lagrange equations. Links with the latter can only be established after 
having decided what are the commutation relations between the dynamical variables, and 
which representation to choose. 

So, which algebra should be satisfied by the dynamical variables? 

Finkelstein suggests {2], 5] that one should quantize in terms of irreducible representations 
of those Lie algebras whose physically meaningful representations are finite dimensional. In 
these notes I want to concentrate on an approach which is, in a sense, complementary to the 
philosophy of Finkelstein (I prefer reducible but infinitely-dimensional representations), and 
whose main ideas can be found already in the two preliminary papers of mine 6|,|7|. However, 
copyright for the term "regularization by quantization" is due to David Finkelstein. 

Restricting the discussion to spinless fields I propose to consider the following two options: 



3 



Canonical commutation relations (CCR), 



[a(p),a( P y) = 5(p,p'), 



(1) 



and harmonic oscillator Lie algebra (HOLA), whose nontrivial commutators read 



[a(p),n{p')} 
[a(p)\n(p>)\ 



8(p,p')a(p), 
-S(p,p')a(p)\ 



(2) 
(3) 
(4) 



with 5(p,p') playing the role of structure constants. 

The idea of replacing Poisson brackets by commutators leads to ([1]) as a natural candidate. 
An interpretation of fields as ensembles of harmonic oscillators suggests (T5])-(j4]). The usual 
identifications n(p) = a(p)^a(p) and I(p) = 1 may seem to imply that CCR determines the 
form of HOLA, but this is not the case, as we shall see in a moment. Moreover, whatever 
algebra one selects, one yet has to tell something about 8(p,p'), which is less obvious than 
one might expect. 

II. CCR DOES NOT DETERMINE HOLA 

Let us now show that a concrete representation of CCR does not yet fix the form of an 
associated HOLA. Let [«i,aj] = l\ be the representation of CCR in a Hilbert space "Hi (l\ 
is the identity map in Hi). Defining rii = a\a± we obtain a representation of HOLA: 



[ai,4] 



I 



(5) 
(6) 



[ai,ni] 





(7) 



Now take 




(9) 



(8) 



One checks that 



[a 2 , 4] = 7i <S>h = h, 



(10) 
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i.e. this is a representation of CCR in H 2 — Hi <8> Hi- In order to extend it to HOLA we 
have to find an appropriate n 2 . Obviously, one choice is 

n 2 = a\a 2 (11) 
= - (a\ai (g 7i + Ii <g a\a\ + a\ (g a x + ai <E> 4) . 

(12) 

The Hamiltonian i7 2 = ftwn 2 represents a system of two interacting oscillators of frequency 
cj/2, where 

H 2 ' = — ajai <g ii + ii (g) — ajai (13) 

is the free part. The interaction term 

H 2 ] = (4 ® ai + a! (g a}) (14) 

is responsible for energy exchange. The algebra is as it should be 

[a 2 ,4] = J 2 , (15) 
[a 2 ,n 2 ] = a 2 , (16) 
[a\,n 2 ] = -a\. (17) 

The problem is that n 2 could be replaced by 

n 2 = n x <g h + h ® (18) 

still satisfying the same algebra 

[a 2 ,a 2 ] = I 2 , (19) 

[a 2 ,n 2 ] = a 2 , (20) 

[4,n 2 ] = -4, (21) 

with the same a 2 , 4 and J 2 . Physically the new representation corresponds to two non- 
interacting oscillators of frequency u if we employ the same correspondence between Hamil- 
tonian and the number-of-excitations operator: 

H 2 = fiwn 2 (22) 
= tuva\ai <g I\ + Ii <g t\jjja\a-i. (23) 



It is not surprising that still further possibilities exist, 



n 2 = a[ ® ai + a\ ® a}, 



(24) 



or an arbitrary convex combination p\n 2 + p 2 n 2 , Pi + P2 = 1- An example of the latter has 
already been encountered, 



We will later see that the above ambiguities imply certain level of freedom in definitions of 
quantum-field 4-momenta and number operators. 

III. NUMBER-OF-PARTICLES REPRESENTATION (BOSONIC FOCK SPACE) 

Before we discuss field quantization in more detail let me introduce here the familiar 
number-of-particles representation of multi particle systems. The construction is interesting 
in itself and is often regarded as the representation appropriate for quantum fields. We will 
later see that the latter statement is not necessarily true, or at least is less obvious than 
what one is generally led to believe. 

One often reads that "there is no quantum field theory, there is only a theory of multi- 
particle systems" . In the next chapter I will give arguments why formal similarities between 
multi-particle systems and quantum fields may be misleading, a subtlety closely related to 
the distinction between particles and quasi-particles. But first we have to understand the 
standard construction, although in the end it will not be employed in my own approach to 
field quantization. 

We begin with a Hilbert space "H whose countable basis will be denoted by \j), j = 1,2.... 
Hilbert spaces that possess a countable basis are called separable (the standard Hilbert space 
from the first semester of quantum mechanics is separable, an example of the countable basis 
being provided by Hermite polynomials). Any operator A has matrix elements (k\A\l) = Am 
and can be written as 



1 1~ 
n 2 = 77^2 + -n 2 . 



(25) 




(26) 



k 



ki 



ki 



If j — 1) 2 . . . is another basis in "H, then 




(27) 



3 
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is unitary, 

= e isoi E = E iSoi*x*i = E = E 1M1 = ^ ( 28 ) 

v+y = El*X*lE ©01 = E l*X*B>01 = El fc >^*01 = El^l = i. ( 29 ) 

relates the two bases, 

v\k) = e©oi*> = Ei^* = i*>' ^ 3 °) 

and has matrix elements 

(*ivii> = (*i E isoio = E<*®^ = ^ 31 ) 

Moreover, 

i*> = m = Ei^>ow> = E^> ( 32 ) 

shows that numbers Vik allow us to write the new basis as a linear combination of the old 
one. 

Now let us take a countable collection of operators Oj satisfying CCR [ajfe,aj] = 5m, and 
let |0) be their common "vacuum state", 

a 3 \0) = 0, j = l,2,... (33) 
(0|0) = 1. (34) 

We do not assume that |0) belongs to H. It is a completely independent object, belonging 
to some linear space J- ', say. Let \j) = a]-|0), (j\ = (0|a 3 -. The new vectors are orthonormal 

(k\l) = (0|a fc oj|0) = (0|a fc aj - aja fe |0) = <5 W (0|0) = <5 W , (35) 

and (by definition) also belong to T . \j) are orthogonal to |0), 

(j|0) = (0|oy|0) = 0. (36) 

Taking numbers V^i = (k\V\l), defined in f[3Tj) . we can define a new set of vectors 

T,m\v\k) = E^w = E^s t i°) = 4io) = i^), (37) 

st 



= (j|V|fc> = (38) 



Note that 

hk = ( E V i^i) ^ = E ^~ ka i = E V kj a ii ( 39 ) 

3 3 3 

so that 

[a*, 4j = [ E v kW E = E 4 E tw°* af] = E v ^ji 

3 l 3 I jl 

= E V kj V im = (V^V) km = (k\V*V\m) = (k\m) = 5 km . (40) 

3 

The fact that 

[a k ,a) m ] = [a fc ,ajj = 5 km (41) 
suggests that there exists a unitary transformation U satisfying 

Ua\tf = 4 = E^a]. (42) 

3 

I will construct U explicitly and show that 

U\0) = tf\0) = |0), (43) 
U\k) = Ua\\0) = UalU^O) = 4|0) = \k) = Y^V jh a]\0) = Y^V jk \j). (44) 

3 3 

As we can see, U will play in J 7 a similar role to that of V in H, 

\k) = U\k) = J2V ]k \j), (45) 

3 

\k) = V\k) = ^2v jk \j). (46) 

3 

Let us now perform the construction of U . To do so we first note that any unitary transfor- 
mation can be written in an exponential way. Indeed, eigenvalues A of a unitary operator 
V, say, satisfy 

(47) 
(48) 
(49) 

A(A|A) = (A|\/- 1 |A) = A- 1 (A|A), (50) 

(51) 



V\X) 


= A|A), 




= A^IA), 


(\\V\X) 


= A(A A), 


(\\V\X) 


= <A|V+|A> 


A 


= A" 1 . 
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Therefore, A = e %tp ^ for some real an unique <p(X) € [0, 2%). Let $ = <^(A)|A)(A|. Then, 

e » = ^ e ^)|A)(A| = ^A|A)(A| (52) 

A A 

commutes with V and has the same eigenvectors and eigenvalues as V. Accordingly, V = e** 
with $"!" = $. 

In particular, there exists X such that U = e x , W = e~ x . It is known that for any X 
and Y 

e x Ye -x _ y + [X,Y] + [X, Y]] + ±[X, [X, [X, y]]] + . . . (53) 

(Of course, provided products and commutators of operators are meaningful, which is not 
obvious if the operators are unbounded.) In order to prove this, let us consider the family 
of operators 

Y(t) = e tx Ye~ tx } (54) 
y(0) = y (55) 



Differentiating (1541) we find 

d 



Y{t) = Xe tx Ye~ tx - e tx Ye~ tx X = [X, Y{t)]. (56) 



On the other hand 



j t {y + t[X, Y] + £[X t [X, Y]) + |[X, [X, [X, Y))} 



t 2 



[x, y\ + t[x, [x, y]] + -[x, [x, [x, y]]] + • • • = [x, Y(t)\. (57) 



We conclude that 

d :Y(t) = j t (Y + t[X, Y] + t [X, [X, y]] + |[X, [X, [X, Y}]] + ...), (58) 



where the differentiated objects satisfy the same initial condition at t = 0. By uniqueness 
of solution of a first order linear differential equation one arrives at 

Y(t) = e tx Ye~ tx = Y + t[X, Y] + |[X, [X, y]] + t [X, [X, [X, y]]] + . . . (59) 

for any t. Putting t — 1 we end the proof. 
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Now take X of the form 

X = ^2 a l x ki a i- (60) 



ki 

Then 



k 

[X, [X, QijW ^ ^ ^ y \ o\i%nrnQ'm-i Qfc-^fcj] ^ ^ ^nm^mfe^fcj ^ ^ a n X nm%rnj 

k nm k nm nm 

= E4(*% ( 62 ) 



Repeating the procedure, we find 

Ua]rf = a] + [X, a}} + ±[X, [X, a}]] + [X, [X, a}]]] + . . . 

= E a * + ^ + ^|( a;2 )fci + ^(^ 3 )fei + • • • ) 

k 

= E4(e%- (63) 

It follows that in order to construct U we first have to express V as V = e x for some 
antiunitary x = —x\ The matrix elements Xkj = (k\x\j) are then employed in construction 
of X satisfying U = e x . Since X|0) = we get 

U\0) = (l + X + ix 2 + ix 3 + . . . ) |0) = I|0) = |0). (64) 

In consequence, once we have a,j corresponding to some basis \j), we know how to construct 
dj such that \j) is related to \j) in a given way, determined by the link between \j) and 

The vectors \j) in J 7 so constructed are completely analogous to \j) in H. Having any 
vector 

we can construct 

ivo = I>w = E^°3 1°) = *) G ^ ( 66 ) 

whose properties are "identical" to those of \ip). ip is an example of field operator. 
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It is interesting to check the action of X given by ( 160]) on i/j\0), 

Xip\0) = ^alxkiai^ipja^O) = J^a^ViK, at] |0) = ^4 x fcrt^<y0) 

fci j klj klj 

= J2 a l x krti\0) = ^0W0*4l°) = 5|0). (67) 

fcZ k 

We find X^|0) = x^|0) but this is true only in action on |0). The general rule for X and ij) 
of the above form is 

[X,i>] = xlp. (68) 

If ip\ is an eigenvector of x, i.e. xip\ = then the .F-space counterpart of the eigenvalue 
problem is 

[XM = AV^a, (69) 

implying 

[x,4> Xi ...$ Xk ] = (\i + --- + \ K )4>\i---4>\ K , (70) 
X$ Xi ...$ Xk \0) = (\i + --- + x K )i>x 1 ...i>x K \o). (Tl) 

The latter formula illustrates the essential difference between T-L and J 7 , namely the fact that 
in T it makes sense to consider vectors of the form 

\jk) = at4|0), (72) 
\jkl) = a|4a||0), (73) 



all of them orthogonal to \l), I = 



acts on these vectors as follows 



0,1,2,.... The operator 



n|0) 


= o, 


n\j) 


= 


n\jk) 


= 2\jk), 


n\jkl) 


= 3\jkt), 
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(74) 



(75) 
(76) 
(77) 
(78) 



so it simply counts how many times one acted with creation operators on |0). For this 
reason n is termed the number-of-particles operator. The subspace consisting of vectors 
satisfying n\ip) = is called the one-particle sector of the (bosonic) Fock space J 7 . This is 
the subspace of J 7 that is "identical" to the initial Hilbert space H we have started with. J 7 
is a countable direct sum of n-particle sectors, each of which is in itself a separable Hilbert 
space. Accordingly, J 7 is also a separable Hilbert space that could become a starting point 
for yet another Fock space — and so on and so forth. 

In 1925 Heisenberg and his coworkers did not know the notion of a Fock space when they 
wrote the first paper on field quantization. In fact, they did not even know quantum me- 
chanics. They did not hear of Hilbert spaces or Schrodinger's equation. But they formulated 
cornerstones of modern quantum field theory, with all its advantages and defects. 

Digression: It is instructive to consider also the following construction. Let us denote by 
% n the Hilbert space spanned by tensor products 

\jl,---,jn) = \ jn), (79) 

and let A n and S n be the projectors on subspaces of, respectively, anti-symmetric and 
symmetric states in % n (for n = we take Ho = C; A = So and A\ = S± are identity 
operators in TL and "Hi, respectively). Let T-L = 0^ o "H n and <x,-,aj : H — > H be defined 
by 

<Xj\jl,---,jn) = C n+1 \jji,...,j n ), (80) 

®j\jo,jl,---,3n) = c n+1 5 jjo \j 1 ,...,j n ), (81) 

aj\0) = 0, (82) 

at|0) = c ± \j), (83) 

= Ci^jO). (84) 

c n+1 and c n+ i are constants to be specified later. Let A = ®^L A n , S = S n , and 

Oj = SaijS, (85) 

a] = Sa]S, (86) 
bj = AajA, (87) 

b) = Aa)A. (88) 
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Let us check: 

aj a]\0) = aj\j) = ci|0), 
a] aj \0) = 0, 

[%,4]|0) = ci|0>, 

1 11 

<¥4b'i) = ^^-(Ij, ji) + j)) = czSaj-djj!) + \ji,j)) = c 2 c 2 -{\ji) + ^Ji)) 

a)a,j\h) = c 1 5 jh a^Q)=c 1 c 1 8 j j 1 \j), 

\ a ^ a ]]\3i) = ^pb'i) + ^Y 5 3h\j) - cici^Jj), 

%a]|ji,j 2 ) = %a^(b'i,i2) + |j2,ii)) = C3a^i(|j)|ji,j 2 ) + |j)|j 2 ,Ji)) 

= c 3 a~(|j, ji, j 2 ) + \j,j 2 Ji) + \j1JJ2) + \j1J2J) + \j2J1J) + \j2,j,jl) 



= C 3 C 3 ^(|ji, j 2 ) + \j 2 ,jl) + S jjl \j,j 2 ) + S jh \j2,j) + S jj2 \jl,j) + Sjj^jji)^ 

= c 3 c 3 ^(\j 1 ,j 2 ) + \ j2J1) + S jjl \j,j 2 ) + Sjji\32,j) + 8m\hJ) + Sjj2\jJi)) 
a)aj\h,h) = a)<Xj\{\juh) + \h,h)) = c 2 Sa]S^(5 jjl \j 2 ) + ^ 2 |ji>) 
= c 2 c 2 S^ (5 jh \j, j 2 ) + S jj2 \j, ji)) 
= c 2 c 2 ^{5 jh I j, j 2 ) + S jh \j 2 ,j) + S jh I j, ji) + 5 jh j>) 

a]] J2> = C3C3- 32) + \ jl) + I j, J2) + ^ | j 2 , j) + <% 2 I jl, j) + <% 2 I j, jl)) 

-C2P2^feib', J2) + 8 jh \j2,j) + <^ 2 |j, jl) + *jj 2 b"l,J>) 

The formulas get incredibly simplified if we take c n = c n = y/n. Then 

[ % ,a>) = |0>, 

\ a ii a ]]\ji) = \h), 

[aj,a]]\ji,h) = ^{\ji,j2) + \j2, ji)) = S\ji, j 2 ). 
Note that since S is a projector then [aj, at] = [aj, Oj]S, hence 

[aj,aj]S\j 1 ,j 2 ) = S\ji,j 2 )- 
Now take any n > 1 and let {cr = (cr(l), . . . , cr(ra))} be the set of all the permutations of the 
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sequence (1, . . . , n). Then 

[% , at] in) = dj a ) ' • • • ' M«) )" fl oE a J l^ 1 ) ' • • • ' M n > ) 



n! 14 — ' J n! 

cr cr 



cr 



n! 14 — ' J n! 

cr cr 



gaj "^T 1 5^ ^l^^ 1 )' • • -'W ~ 5a l"~r X^^'MD ^l^ 2 )' • • 

Vn+1 1 
n! n + 



= Sa 3 '"'^n+l Yl {\j>MV> • • • '^(«)) + • • • + iMl), • • • Ja(n),j)) 



cr 



= S (iM 1 )' • • •>>(»)) + *iJ,(l)b'iM2)i • • • iJ»W> + • • • + ^ CT (i)l>(2), • • • , jcr(n), j)) 

^^Ml)b'> Jcr(2), • • • , jcr(n)) 

a 

= (l^ 1 )' • • -'W + ^Ja(l)b''^(2)> • • • , jcr(n)) H h ^J CT(1) |jcr(2), • • • , jcr(n), j)) 

CT 

~ J2 6 3MD {ij'MVi ■ ■ ■ >Mn)) + • • • + |jcr(2), • • • , Jcr(n), J>) 

cr 

= 5 ~| l^lMl),---,Mn)> 
= S\j!,...,j n ) 

Finally, 

[%,a}] = S (89) 
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Analogously 

bjb]\0) = ^|j) = Cl |0), 

b]bj\0) = 0, 

{b„b]}]\0) = Cl ]0), 

b j b }\ji) = C2bj-(\j,ji) - = c 2 Aa j -(\j,j 1 ) - \j u j)) = c 2 c 2 -(\ji) - 5 jh \j)) 



b ) b i\ji) = c 1 6 jjl b]\0) = c 1 c 1 6 jjl \j), 
{bj,b]}\ji) = ^Y'\ji)-^Y'S jjl \j) + c 1 c 1 6 jjl \j), 

bjb]\ji,j 2 ) = bja^djuh) - |j 2 , Ji>) = CabjA^djJ^jz) ~ \j,h,ji)) 

= c 3 Aa j -(\j,j 1 ,j 2 ) - \jJ2J1) - \j1JJ2) + \j1J2J) - \h,ji,j) + \j2, j, ji>) 

= c 3 c 3 A^(\j u j 2 ) - \j 2 , j x ) - 5 jh \j,j 2 ) + 6 jh \j 2 ,j) - S jj2 \j u j) + <% 2 |j, ji>) 

= c 3 c 3 ^ (\ji,j 2 ) - \j 2 ,ji) - S jh \j,j 2 ) + S jh \j 2 ,j) - S jja \j u j) + SjjAjJ^ 

bjbj\j u j2) = ^-(Ijuh) - U' 2 , Ji>) = c 2 Aa)A-(8 jh \j 2 ) - 5 ji2 |ji)) 

= c 2 c 2 A^(5 jh \j,j 2 ) -5jh\3,h)) 

= c 2 c 2 ^ {5 jh I j, j 2 ) - 5 jh |j 2 , j) - S jh I j, ji) + 5 jh |j'i,j>) 

{bj,b\}\ji,j 2 ) = c 3 c 3 -(|ji, j 2 ) - \j2J1) ~ Sjn\j,j2) + S jh \j 2 ,j) - S jj2 \j u j) + 5 jh ./-./:)) 

+c 2 c 2 ^(8 j:jl \j,j 2 ) - 8 jh \j 2 ,j) - 5 ji2 |j, ji) + <^ 2 b'i, j)) 

As before, choosing c n = c n = y/n, we simplify the formulas. Let Yl' a = So- s ig n (°")- 

1 x 1 . 

{bj , &J } I jl , . . . , j n ) = bj — >a ) , • • • , j<r(n) ) + b) -y ^ |j CT (l) , • • • , Ja(n) ) 



n! 14 — ' J J n! 



/(/. /i. 



= ^ w+1 Z] / (l J ''^( 1 )' • • • ' 3 ^ + ■■■ + (-l) n bV(l), • • • J*(n),j)) 

cr 
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= A ^J2'i} 3a{lh ' ' ' ' 3<T{n) ^ + ( -1 ) 1<J JMl)b''^(2), • • • ,Ja(n)) + ■■■ + (-l) n ^i CT(1) |j CT (2), • • -,ja(n),j)) 
a 

+A ~\ ' • • • ' Mn)) 

= A ~^\ • • "'^W) + ( -1 ) 1<J Wa(i)U'^(2)> • • • .^n)> + • • • + (-l) n ^ ( i)|j'<T(2), • • • ,j*(n), j)) 

S'^'-d) (b''^(2)> • • • . jV(n)) + • • • + (-l) n_1 ba(2), • • -,ja(n),j)) 

cr 

= A ~|l^'l^(i)'---'^W) 

cr 

= A|ji,...,j n ) 

Finally, 

{M5> = A (90) 

In the subspace of symmetric states 5 CH, 55 = 5, the commutator [%,at] = S acts as 
the identity operator, 

[a j7 a]]S = S. (91) 

Similarly, in the subspace of antisymmetric states A C H, AA = A, the anti-commutator 
{bj, bj} = A acts as the identity operator, 

{b J ,b]}A = A. (92) 

The two algebras (canonical commutation and anti-commutation relations) in physical appli- 
cations will correspond to different-spin representations of SU(2) or SL(2,C) groups. There- 
fore formulas such as aj|0) = \j) and &]|0) = \j) should not both generate the same The 
simplest solution is to start with |0,0) = |0) ® |0), 

a]\0) <g> |0) = \j) <g> |0) e S <g> A, (93) 
&}|0)<g>|0) = \0)®\j) eS®A, (94) 



and so on. The algebras then read 



[ aj ,a]] = S®I, (95) 
{bjM} = I® A. (96) 
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This type of construction can be found in the Berezin book [8j but, apparently, goes back 
to the original works of V. Fock. It is interesting that the right-hand sides of the algebras 
are not given by the identity operator but by projectors on subspaces of totally symmetric 
or anti-symmetric states. A 



IV. ENSEMBLES OF INDEFINITE-FREQUENCY OSCILLATORS 

Energies of classical free fields look in Fourier space analogously to those of ensembles 
of oscillators. This observation was at the roots of the first approach to field quantization, 
formulated already in 1925 by Heisenberg, Born and Jordan [9]. One should bear in mind 
that the first paper of Schrodinger, explaining the role of eigenvalues of operators, appeared 



a few months later |l0|. The authors of {9] apparently die 



not yet understand the idea of 



quantum superposition (as suggested by Max Jammer in 11(). Having to consider oscillators 
with different frequencies they basically had no other option but at least one oscillator per 
frequency. In this respect nothing essential has changed in mainstream quantum field theory 
since 1925. 

I will now show that their basic assumption is not at all natural. To do so, let us consider a 
simple 2D pendulum in linear approximation. Classically, its ground state would correspond 
to no oscillation at all. Quantum mechanically it would imply vanishing momentum and 
fixed position, a possibility excluded by the uncertainty principle. In consequence one finds 
the ground state oscillation with energy \ fi\fgjl = \tvjj{l) (I is the pendulum length). It also 
practically means that the "lowest atom" of the pendulum is described by the center-of-mass 
wave packet ^(X) where |\1/(X)| 2 is a Gaussian. 

If the pendulum is suspended at the origin (X, Y) = (0,0), the length is given by / = 
\Y\ = —Y. For a true pendulum the 2D wave function \I/(X, Y) would have to be smeared 
out also in the Y direction and, in fact, a more realistic model should employ a nonseparable 
potential U(X,Y) = mu(Y) 2 X 2 /2, with deep conditional minimum at Y — —I. Note that 
the frequency is no longer a parameter but an eigenvalue uj(Y), in position representation, 
of some operator Q. Similarly to Schrodinger's cat that exists in superposition of being dead 
and alive, our pendulum exists in a superposition of all its possible lengths, so that many 
different frequencies can be associated with a single oscillator. The example is generic - 
in sufficiently realistic cases the us are not just classical parameters but functions of other 
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observables. 

I don't see any reason why at levels as fundamental as those related to field quantization 
the appropriate us and wave vectors k should be more classical than Q from the preceding 
example. I would rather expect the wave vectors and frequencies to be eigenvalues of some 
operators. In order to understand formal implications of the latter postulate we first have 
to understand the simplest examples based on nonrelativistic oscillators. 



A. Single indefinite-frequency oscillator 

The special case Q = ul, with parameter oj and identity operator I, is equivalent to the 
standard harmonic oscillator. 

The simplest nontrivial generalization of Q, — oul occurs if the operator Q has discrete 
spectrum and commutes with canonical momentum and position. So let p and x be the 
canonical momentum and position acting in the Hilbert space spanned by the number-of- 
excitations eigenvectors \n). Now consider the following representation 

ft = ^2<jj\w)(u}\®I, (97) 
P = ^\uj){uj\®p = I n ®p, (98) 

Ul 

Q = ^2\u)(u\®x = I n ®x, (99) 

U! 

with Hamiltonian of the usual form 

P 2 mn 2 Q 2 Ml ( + t \ / t 1 \ . . 

2^ + ^^ = ^r\ aQan + ann ) = hn { a n a n + 2 In ® 1 )' (100) 

and CCR 

a n = 1 (mQQ + iP), (101) 
\/2hmVL 

at = - 1 (mttQ-iP), (102) 

acting in the Hilbert space spanned by \u>,n) = \u) <g> \n). For any u 

H\u,n) = tvjj(^n-\ — j|w,n), (103) 
a n |w,0) = 0. (104) 
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Let 



1^) = ^2*Pu>,n\uJ,n} 

u),n 

be an arbitrary state. The average 

(^\H\^) = J2 nw { n+ l)^ 



(105) 



(106) 



looks as an average energy of an ensemble of harmonic oscillators, different frequencies 
occurring with probabilities 



Pu> / , I ' ,;.n 
n=0 



We can also write 



oo 1 



U! n=0 



Here 



\/2hmuj 
1 



v 7 ^ 



mwi + ip) 
(mcux — ip) 



mw 



-x + 



2£ v 7 ^ 



=P, 



muj 



-x — 



2h 



--p, 



x = 



2h 



a w + a 



are the usual creation and annihilation operators, and 



a ul = 



\u){u\ ® a u , 
\u)(u\ ® a*, 



V V / 



(107) 



(108) 



(109) 
(110) 
(111) 



(112) 
(113) 
(114) 

(115) 
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We have arrived at the decomposition 

H = J2 H ^ ( 116 ) 

H w = ^(ala^ + a^al) (117) 

= hujn^ + ^-I^ (118) 



oo 



(1>\H U \1>) = hwJ2(n+^)\^ n \ 2 . (119) 

n=0 

As we can see, H w has properties of Hamiltonian of an oscillator whose frequency is u, and 
if is a sum of such Hamiltonians taken over all the eigenvalues of f2. The new feature is the 
fact that H describes a single harmonic oscillator existing in superposition of different u>s. 
In Heisenberg picture we find 

a u (t) = e im l h a u e- iHt ' h = (120) 
a w (ty = e iHt/R ale- im/h = a* e*", (121) 



m = V^(^ e "" + a - e ")' (122) 

Q u (t) = V^( a - e "^ + °« e<wt ) ' (123) 

which are again the standard formulas. However, in spite of all these similarities to the 
known textbook results, the representation of HOLA is here a nonstandard one: 

= alau = \u))(uj\ <g> ala u , (124) 

I u = \u){u\®I, (125) 

[a^alA = Sufi/Iu, (126) 

[a w ,rv] = 8 Ui0J >a u , (127) 

l a i, n u>>] = -<W' a Ir (128) 

I u commutes with all the elements of the Lie algebra, but is not proportional to the identity. 
By Schur's lemma the representation is thus reducible. Its irreducible representation com- 
ponents are spanned by \u>, n), with fixed u> and arbitrary n. The projector on this subspace 
is given by I u itself. 

Physically, the above representation describes a single oscillator whose states are wave 
packets consisting of superpositions of different eigenvalues u>. 
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Remark: Let us note that (11091) and (11101) imply 



-. (mux + ip) , -^=^= (mu'x — ip) 

y2hmu y2hmu' 



= [mux, —ip] H [ip, mu'x] 

2hm\/uu' 2hm\/uu' 

= ^±£ = L (u,u'). (129) 
2yuu 

The right side of fl 1 2 9 [) is a function satisfying 

l(u,u') = l(u',u), (130) 
i(u,u) = 1, (131) 

which is neither 1 nor any kind of delta. The explicit form of l(u, u') does not occur in HI 26[) 
due to the presence of in HI 12[) and (11131) . We should keep this observation in mind 

when we perform field quantization. A 

B. Several indefinite-frequency oscillators 

Now, what about two such oscillators? One can trivially extend all the operators by 

->• a u ®I(l), (132) 

a, ->• 1(1) ®a u , (133) 

and so on, where 

Hi) = J2^ ( 134 ) 

is the identity in the one-oscillator Hilbert space spanned by \u,n). However, if one ad- 
ditionally requires their bosonic statistics, something one expects for spin-0 fields, say, the 
operators should preserve symmetry of states. The natural bosonic generalization is then 

a u (2) = c 2 (a u ® 7(1) + 7(1) ® aj) (135) 

where c 2 is a constant, and 



[a w (2), a< y(2)t] = ^|c 2 | 2 ^®7(l) + 7(l)®7 u 

= S U ^I U (2). (136) 
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Let us note that 

£}/ w (2) = 2|c 2 | 2 7(l)®7(l) =2|c 2 | 2 7(2) (137) 

suggesting the normalization c 2 = l/y/2, analogous to (JSj) . In order to get the whole HOLA 
we have to define n,(2). The first guess 

n,(2) = ^(2)^,(2) (138) 

leads to 

[a,(2),n,(2)] = <L,u/7,(2)a,(2) (139) 

which is not even a Lie algebra. (Certain observables do satisfy such generalized Lie algebras 
- Hamiltonian, angular momentum, and the Runge-Lentz vector of the Coulomb problem 
provide an example.) But we remember that in standard representations of CCR we have 
encountered the ambiguity of ( ITT]) versus ffTHj) . In the present context the second option 
reads 



n,(2) = n w <g> 7(1) + 7(1) <8> n„ (140) 



leading to HOLA 



^(2)^^(2)+] = 5 U ^I U (2), (141) 

[a w (2),7v(2)] = 5,^0,(2), (142) 

[o w (2)t,n u ,(2)] = -<J WlW ,o w (2)t. (143) 

This representation is again reducible since 

7,(2) = i(/ w ®J(l) + J(l)®J w )^/(2), (144) 

J>(2) = J(2). (145) 

71,(2) counts the number of excitations of the two-oscillator system. Especially interesting 

is the form of 7,(2): This is the number-of-successes (in two trials) operator known from 

works on probabilistic interpretation of quantum mechanics 12Hl6|. Eigenvalues of 7,(2) 



are 0, 1/2, 1 and describe the fraction of positive answers (in two trials) to the question "is 
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the frequency of the oscillator equal to uT\ if u is selected at random. Due to the resolution 
of identity f!145j) the element 7 W (2) is a positive operator- valued measure 

An extension to arbitrary number N of oscillators is now clear, with = 1/yN. The 
fact that Iu(N) becomes the iV-trial number-of-successes operator is essential for the limit 
N — > oo which can (and later will) be treated by weak laws of large numbers. 

Before I make a digression on scalar fields let me note that in Heisenberg picture, for 
arbitrary N > 1, we find 

a u (t,N) = e^^ n a u (N)e^^ A = a u (N)e- iut , (146) 

if the Hamiltonian is given by H(N) = ^ w hio h u (N). For H(N) = fkv a u (N)^ a^N) we 
would get 

a„(t,N) = e iH{N)t/h auJ (N)e- lH{N)t/h = a w (N)e~ iulI ^ N)t . (147) 

This is a strong argument in favor of H(N) in contrast to H(N)... (Still, as mentioned 
already in I have doubts here — the choice of H(N) is appealing for various reasons, so 
we will devote some time also to this option.) 
Canonical position operators 

QUt,N) = .[^(a^e-^ + a^NYe 1 ^), (148) 
V mu V / 

Q(t,N) = J2J^(aUN)e-^ + a u (Ny e ^y (149) 

LU 

are formally very similar to scalar-field operators. This observation is behind the basic 
physical intuition that leads to quantum fields, quantized in "my" way. 

V. FREQUENCY-OF-SUCCESSES OPERATOR AND HOLA 

For N = 1 the operators 11^ = Iw an d Tlffl = 7(1) — Id are orthogonal projectors : 
(nL 0) ) 2 = nL 0) , (U^) 2 = n^, nL 0) nL 1} = 0, n£ 0) + nL 1} = 1(1). The central element 

I U {N) = i (l u ® 7(1) ®---® 7(1) + + ®7(1)® /„) (150) 
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has eigenvalues 0/N, 1/N, 2/JV,..., (N — 1)/N, N/N. Its spectral decomposition can be 
deduced from 

1 



Hi) = nl 0) + nw = ^nw (151) 

s=0 

I U (N) = ^ (nW <g> 7(1) <g) ... <g) 7(1) + • • • + 7(1) <8> • • • <g> 7(1) ® n^) 

= ^n(°)®...®n(°) 

+ ^ (n« ® 7(1) <g> • • • <8> 7(1) + • • • + 7(1) <g) ... (g) 7(1) (g) n«) (152) 



and 



Since 



7 w (7V)n^) g) ... (g) n£") = Si + " n + Sn u^ ® • • • <8) n<>">. (153) 



1 

5^ njw <g> • • • <g> n[ sjv) = /(iv) (154) 

si,...,sjv=0 



we arrive at 

1 



UN) = J2 UN)n { : i] ® • • • ® n<?"> 

si,...,sjv=0 



1 TV 
Si H h Sat, 



E ^^V ^..^)^^), (155) 

si,...,sjv=0 s=0 

n -(^) = E n&° <s» • • • o n&^. (156) 

siH hsjv=s 

Now let 1) = E W) n^«,nk,n), N) = 1) <g> • • • <8> 1). The average 

(V, 117^, 1) = linL 1 ^, 1) = 5^ l^.nl 2 = P«, (157) 

n 

is the probability of finding u. The average 

siH hsjv=s 

= £ (^,i|n^|v,i)...(v,i|n(;-)|v,i) 



SlH hSAT=S 

/ N^ 
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is the probability of finding u exactly s times in N measurements, performed on each of the 
iV oscillators once, if each oscillator is in state \ip, 1). 

If F : [0, 1] — > R is continuous, then by the weak law of large numbers for binomial 
distribution (or, more generally, Feller's theorem jl8|) 



Km {1/>,N\F(I U (N))\1/>,N) = lim (V, N\ V F (|- ) H 

N^too TV— »oo z — ' \iV/ 

s=0 



N 



iV) 



N 



lim X^f(— 

N^oo^ \N/ \ S 
s=0 



N 



(159) 



So, in practice, predictions concerning the gas consisting of iV indefinite-frequency oscillators 
are in the limit iV — > oo close to those of the standard oscillators with the modified CCR 
[ a oj,^,] = Pu5 ww i. This result will be the basis of the correspondence principle relating 
standard regularized quantum field theory with the one I propose in these notes. In fact, as 
discussed in detail in 
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2lj . there exist two physically meaningful limits iV — > oo. One is 



just the weak law of large numbers while the second one is a kind of thermodynamic limit. 



VI. STATES OF INDEFINITE-FREQUENCY OSCILLATORS 

States of indefinite-frequency oscillators have properties that will be later used in con- 
struction of vacuum, n-particle and coherent states of quantum fields. 



A. iV-oscillator analogues of n-photon states 

The "vacuum subspace" consists of vectors that are annihilated by all annihilation oper- 
ators. The iV-oscillator vacuum subspace is spanned by 0) <g> • • • <g) \cun, 0). I will assume 
that vacua are given by pure product states 

\0,N) = Yl O^-OJwi.Ojg-^kO), (160) 

£i»l,...,u;jv 

where the O^s play a role of a single-oscillator wave function, normalized by 

Ei°-i 2 = L ( 161 ) 

Puj = \Ouj\ 2 is the probability that a given oscillator has frequency to. 
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Let us begin with the analog of an ordinary 1-photon state, that is a u (NY\0, N). Its 
squared norm is 

(O^M^iN^N) = (0,N\I u (N)\0,N) 

= (0,1|/ W |0,1) 

= \OJ 2 = PuJ . (162) 

Whatever representation of n^i^N) satisfying HOLA [n u (N), £v(iV)t] = 5 tMJ ia u i(N)^ and 
annihilating \0,N) we take, we find that a u (Ny\0,N) is a "single-excitation" state, i.e. 

MiVK(iV)t|O,A0 = [n u (N),a u (Ny]\0,N)=a u (Ny\0,N). (163) 

The peculiarity of this representation of HOLA is that there exist "single-excitation" states 
that are not spanned by vectors a u (N)^\0, N). A simple (and generic) example is provided 
by any state of the form 

F(UN)) auJ (Ny\0,N). (164) 

In general, n-excitation states are any states spanned by 

|o>i, rii) <g> • • • <g> \u>n, un), ni + h n N = n. (165) 

Their particular sub-class is given by 

a Wl (iV) t ...a Wn (iV)t|0,iV). (166) 

Indeed, each action of a creation operator adds one excitation to the state, and \0,N) has 
zero excitations. 

B. iV-oscillator analogues of coherent states 

Quantum mechanics textbooks give (at least) two definitions of coherent states of a 
harmonic oscillator: States generated by displacement operators acting on the oscillator 
ground state and eigenstates of annihilation operators. In reducible representations of HOLA 
the two definitions are not equivalent. 
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The displacement operator 

= exp ( - i ^ \<*u\ 2 Iu( N i) exp ^ ^ a u a u (N)^ exp ( - ^o^a w (iV)j. 

(167) 

is unitary 

£>(o:,A0 t = Via^y 1 = V(-a,N). (168) 
The operator can be also represented in the following useful way 

V(a, N) = 6^ E " (au-Mljt-Q^Ml)) $ . . . (g, e 7F (^Ml^-c^Ml)) 

= V(a/VN, 1) <8i ... (8) T>(a/VN, 1) 

= V(a/y/N, lf N , (169) 
Z>(a, 1) = exp ^ ^ <g> - 5~a) j 

= ^ <H> exp (a^flt — a^a) 

a) 

= ^ ® exp ( - \a UJ \ 2 /2) exp (a a; a t ) exp ( - a^a) 

= J2\u)("\®V( aiJ ). (170) 

V^a^) is the ordinary displacement operator known from quantum mechanics textbooks. 
The name of the operator comes from the following property 

V(a, Ny auJ (N)V(a, N) = a^N) + aJ w {N), (171) 
V{a,N)^a w {N)^V{a,N) = a^N)^ + c£l u (N), (172) 
V(a,N?I u (N)V(a,N) = I U (N). (173) 

Combining these formulas we obtain a generalized eigenvalue problem 

a w (N)V(a,N)\0,N} = V(a, N)V(a, N) jf a UJ (N) / D(a, N)\0, N) 

= V{a,N)(a u {N) + a w I w {N))\0,N) 
= auV{a,N)I„{N)\0,N) 

= a u I u (N)V(a,N)\0,N). (174) 
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The coherent state 

\a,N) = V(a,N)\0,N) (175) 
= exp(-i^|« w | 2 4(iV))exp(^a^(iV)t)|0,iV) (176) 

LU LU 

= V{a/VN, l)\0, 1) ® • • • <g> V(a/VN, 1)\0, 1) (177) 

= |a/V^V, 1) <S> • • ■ ® |a/V^V, 1) (178) 

satisfies 

auJ (N)\a,N) = a u I u {N)\a,N). (179) 

The latter is a generalized eigenvalue equation in the sense that it combines several ordinary 
eigenvalue problems that can be extracted from (11791) by means of Uu(s/N): 

a w (N)n w (s/N)\a,N) = U w (s/N) auj (N)\a, N) 

= ^aMs/N)\a,N). (180) 

Vector H^s /N)\a>, N) is the eigenvector of a w (N) with the eigenvalue sa u /N. 

The difference with respect to the "ordinary" harmonic oscillator known from textbooks 
is that the latter involves I w = 1 so that the coherent state is just a single eigenvector with 
eigenvalue a u . We will later see that the presence of s/N may make various quantum field 
theoretical predictions more physical (if quantum fields are defined by means of an analogous 
representation of HOLA). 

C. Statistics of excitations 

The n-excitation subspace is spanned by vectors ( I165p . The projector on this subspace 
is given by 



U(n,N) = yj \ui,ni)(u)x,ni\ <g) ■ • • <S> \uj N ,n N )(u) N , n N \ 

mH hrtjv=™ LUi...LUff 

( J ^® KK^il) (in® \n N )(n N \)- (181) 



JtiH hnjv="- 

Probability of finding n excitations in a 1-oscillator coherent state is 



\2\n 

|2 



p(n,l) = (a,l|n(n,l)|a,l) = ^|0^| 2 (|Q; "| ; e ~^M\ (i 8 2) 
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(1182p is the Poisson distribution typical of standard coherent states weighted by probability 
that frequency of the oscillator equals lo. This is precisely the result one expected for a 
single oscillator wave packet. 

Probability of finding n excitations in an iV-oscillator coherent state can be computed if 
one recalls that f)166p belongs to the n-excitation subspace. Therefore, 

p(n,N) = (a,N\U(n,N)\a,N) 

( ^(0,iV|exp(-^|«J%(iV))(^c^a Wl (iV))"(^a W2 a W2 (iV)t) n |0,iV) 



Taking into account 



[y^QimO^^iV), (^a^a^jN)^ = ^ ja^ (N) , ( ^ a W2 (N) ' 



n 



0)2 

(N) 



n-l 



and acting with 



UJl U)2 



Wl 
n-l 



W2 
in 



+ [(^c^MiV))" , (^ a ^ a ^W t ) n ] ^SSW (183) 



on |0,iV}, one finds 

' \ Tl / — — \ TL 



n 



/ \ n— 1 x \ n— 1 

E k i 2 / wi (ao ( ^ a -i ( E a ^ i°' ^ 



n!(£|aUX(w))Vw>. 
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The end result 



p(n,N) = (0,AT|exp(-^|a w |%(iV))i(^| awi | 2 / tJl (iV)) n |0^) 



1 d n 



-(0,N\exp(\J2M 2 UN))\0,N) 



A=-l 



(184) 



is a generalized Poisson distribution. In order to better understand the generalization we 
have found let us take a closer look at the generating function 



(0,N\exv(\J2M 2 UN))\0,N) = (0, 1| exp (a^ -^K| 2 I w ) \0, 1) 

Ul Ul 

= (0,l|^^exp(lA|aJ 2 )|0,l) 



N 



N 



Ei°. 



N 



The possibility of taking the sum in front of the exponent comes from the fact that I w = 
\u>)(oj\ ® 1 is a projector. Now, let us introduce a new parameter q satisfying 1 — q — l/N, 



q)X\a u \- 



2 \ 1-q 



expln^^lOjV 1 -^! 2 ) 1 - 9 

exp ^ I ^ln(^|O w |V 1 -«) A ^l a )j 



(185) 



Expression under the exponent, 



r ^ln(ElO»|V 



-<?)A|a„| : 



186) 



is well known in probability and information theory: This is the Kolmogorov-Nagumo av- 
erage of the random variable Aja^l 2 [2^|, introduced by Alfred Renyi in his derivation of 
generalized entropies [23|. Since 



hmln^lOjV 1 "^! 2 ) = ln(^|OJ 2 ) = 0, 



(187) 



(0,1[0,1)=1 



we can use the de l'Hospital rule to compute the limit 



!r l-in(5>jy'-^) = Ay> 



lim 

i — q 



|2i |2 

a; Ota 



(188) 
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The limiting generating function 

Ul ' \ ^ Ul 

(A^|OJ>J 2 ) (189) 



= exp 

generates the Poisson distribution 

p(n,oo) = ^exp(A^|0.|>.| 2 )| A=i . (190) 

Ul 

Asymptotically, for large N, the gas of coherent-state indefinite-frequency oscillators pos- 
sesses Poissonian statistics of excitations. What is interesting the parameter of the Poisson 
distribution is not just |a: w | 2 but (O^Ha^l 2 , a result of great importance for my formulation 
of quantum field theory. For finite N the distribution is a Renyi-deformed Poissonian. 
Let me summarize this section. For finite N we get 

(O, N\ exp (A M 2 UN)) \0,N). (191) 

Ul 

The N — > oo limiting case is 

exp (a J^p w |a w | 2 ), p u = \O u \ 2 . (192) 

Ul 

The standard formalism based on infinitely many standard oscillators would imply 

exp(A^|a w | 2 ) (193) 

Ul 

which, in practice, is often replaced with 

exp (A^xJcU 2 ), (194) 

Ul 

where < Xw < 1, hm w ->oo Xw = 0, is a cut-off function introduced by hand if |« w | 2 = oo. 
So our has automatically appeared in the place where standard formalism is artificially 
amended by adding Xu>- 

Remark: The limit q — > 1 is known in information theory as the Shannon limit because the 
Kolmogorov-Nagumo average of random variable ln(l/|0^| 2 ) (the amount of information 
obtained by observation of an event whose probability is |OJ 2 ), 

^In^lOjV 1 -^!^! 2 )) = _^l n ^|O w |V n l -l 2(<, - 1) ) 

^ Ul ^ Ul 

~in(E(m 2 )') 

Ul 
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tends to Shannon's entropy 

^ iT ^\n(^\O u \ 2 e^^°^) = £|OJ 2 ln(l/|OJ 2 ) = -]T|OJ 2 ln|0 



|2 



Generalized entropy -A- In (^Y^^p'Cj * s termed the Renyi g-entropy of probability distribu- 
tion p w . Similarly to Shannon's entropy it is additive for independent events 



= t^(£<) + t^>"(£^)- 



q \^ V l-g 



A less obvious property of the Kolmogorov-Nagumo-Renyi average is the following analog 
of 

(A + C) = (A) + C, (195) 
for a constant C and random variable A: 

(A + C) q = ^^(^p/^^^) 

= r L 7 ln (Ep^"- 4 ') + r L 7 ln K ,,c ) 

= (A), + C. (196) 

General Kolmogorov-Nagumo averages corresponding to a monotonic function are defined 

as 

( A )* = _1 ($>u,<KA,)). (197) 
It is interesting that only for exponential or linear one finds (A + C)^ = (A)^ + C.k 

VII. DIGRESSION ON FREE SCALAR FIELDS 

Consider a free scalar field of mass m. The latter means that in momentum (i.e. Fourier) 
space the field is defined on the manifold of 4-momenta p satisfying the constraint p 2 = p\ — 
p 2 = m 2 (in sections dealing with relativistic fields I work in units where c = 1 and h — 1). 
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Due to the constraint the four components of p are not independent: p = ±a/p 2 + to 2 . 
The manifold is a 3-dimensional hyperboloid in M 4 , consisting of two sheets corresponding 
to the two signs of p . Depending on the sign we speak of hyperboloids of, respectively, 
future-pointing (p > 0) and past-pointing (p < 0) 4-momenta. If m = the hyperboloid 
is termed the light cone. 

Lorentz transformations are linear transformations p' = Lp that do not change p 2 ^ 0. 
If p 2 = 0, with p 7^ 0, the (conformal) group of transformations is larger than the Lorentz 
group. It contains, in particular, rescalings of the form p' = 5p and the like. If p = then all 
linear transformations preserve p 2 . For any Lorentz transformation detL = ±1, a condition 
implying that dp' d 3 p' = d 4 p' = ±d 4 p = ±dpod 3 p. 

Now assume that p = a/p 2 + m 2 > 0. Changing variables from (po,p) to {m,p) we find 
that 

d A p = dmd 3 p^^^= = dim 2 ) — P (198) 
^Jp 2 + m 2 2\Jp 2 + m 2 

Since {p') 2 = (Lp) 2 = p 2 = m 2 it follows that 

d*P = d(m 2 ) =±rfV 

= ±d(m 2 ) (199) 
2 a/p' 2 + m 2 



and thus (for any to, even m — 0) 



(200) 



2^p 2 + m 2 2v/p' 2 + m 2 ' 

if the /our components of p and p' are related by a Lorentz transformation. The same 
argument can be applied on the past-pointing part of the mass-m hyperboloid — just change 
variables from (po,p) to (—m,p). 

In consequence, in any integral we can change variables according to the recipe 

d3 " ; F(p) = / ^L^ F(P >) = f ^L=F { L P) . (201) 



2a/p 2 + m 2 Jrs 2a/p /2 + m 2 Jr3 1\fp 



m 



The Lorentz invariant measure on mass-m hyperboloid (divided by (2n) 3 for certain reasons 
related to Fourier transforms) will be denoted by dp, i.e. 

dp = d ] P (202) 

(27r) 3 2vV + m 2 
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In many applications it is convenient to work from the outset with 3D integrals involving 
dp and not just d 3 p. Now let px = PoXq — px (I will sometimes denote px by p ■ x or p^x 11 ). 
Two arbitrary solutions 

<j> x {x) = y"dp(a 1 (p)e-^ + 6 1 (p)t e ^), (203) 
2 (» = J ' dp(a 2 (p)e~ ipx + 6 2 (p) t e* a '), (204) 
of the Klein-Gordon wave equation 

(□ + m 2 )0(a;) = 0, (205) 
lead to a family of conserved Noether currents, <9 M T 'J = 0, 



T/(x) = C^d r Mx)d»Mx) + d^M^WMx)) (206) 



We interpret a(p) and 6(p) as amplitudes describing particles and antiparticles, respectively. 
C is a normalization constant [C = 1 for charged fields, C — 1/2 for neutral fields i.e. those 
whose particles equal to their antiparticles: a(p) = b(p)}. The 4- vector 



J d 3 xT^ (x ,x) (207) 

= j ' dpp fl {b 1 ( P ya 2 (p) + a 1 (p)b 2 (p)^ (208) 
is time- independent. Another conserved current is 

jp{x) = iq<j>i{x)dn<p2(x) ~ iqdp<pi(x)fa(x) (209) 
with the time-independent scalar 

q[<t>i,<h] = J d 3 x J (x ,x) (210) 

= q J dp{b x {p)^a 2 {p) - ai (p)b 2 ( P yy (211) 

The transition from (12071) to (I208p . and from (12101) to (12111) . involves the assumption that 
it is justified to use the Fourier transform representation of Dirac's delta, 

5 (3) (P) = 7^53 J tPxe**, (212) 
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which I will comment on later. (Of course, one could start with (12081) . (121 ip . and not with, 
(I207p . (12 1 Op : then the comment on applicability of (121 2p would be irrelevant.) 
Inserting 

4> 2 {x) = 00)= / dp(a(p)e- ipx + b( P ye ipx ^ (213) 
(j) X (x) = 0(x) f = J dp(b{p)e- ipx + a(p)^e ipA j, (214) 
we find the 4-momentum 

P M [0t,0] = Jd PPft (a(pya(p) + b(p)b(p)^ (215) 

and charge 

qWA] = Q J dp(a(pMp) -b(p)b(p)^. (216) 

What one usually finds in the literature are the particular cases (I213p and (I214p . while the 
general ones are not mentioned. The general expressions are complex (or non-Hermitian, if 
quantized) and this is probably why they may seem "unphysical" , although their real and 
imaginary (Hermitian and ant i- Hermit ian, respectively) parts are separately conserved. 

However, when one comes to field quantization it may pay to maintain their general 
forms. Let me give several applications. 



A. Bosonic Fock space 

Let us begin with the orthodox formulation based on bosonic Fock space representation: 
a(p) = ajipj(p), b(p) = ^2jbjipj{p). cij and bj are bosonic annihilation operators and 
ipjip) are orthogonal functions 

J dp^{p)Mp) = S jk . (217) 
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Then 

P^\<P\ = J d 3 xT^(x ,x) (218) 
= j dpp lx (a(p) j 'a(p) + b(p)b(p)^ (219) 

= E / dpp v (^3(PWk(p)a}a k + 4)k(p)4>j(p)bkbf) (220) 
jk J 

= I d PP» ^(P)Mp)l(a]a k + b k bfj (221) 
jk J 

= J2Wjk(a]a k + b k bty (222) 

jk 

g[0 f ,0] = y rf^Jo^o,^) (224) 
= 9 1 dp(a( P ya(p) - b(p)b(p)^ (225) 

The "vacuum terms" X^-^/Jjj and — 1 are divergent. The infinite negative charge of 
vacuum may be intuitively interpreted as the charge of the "Dirac sea" of antiparticles. The 
remaining charge term 

:q: = q- g vacuum = <?E ( a i a i " h ) h l) = ~ n -) ( 227 ) 

3 

looks reasonable. 

A similar analysis can be performed for the case of 4-momentum. For example, the 
Hamiltonian H = P leads to vacuum energy 

^vacuum = E / dp^p 2 + m 2 \^(p)\ 2 > I dp\^(p)\ 2 = m^l. (228) 

3 3 3 

On the other hand 



H . H ^vacuum 



[a)H jk a k + &tif jfc & fc ) (229) 

jk 

has the form we know from general considerations on the Fock space J 7 . In spite of these 
partly acceptable predictions the presence of infinite terms makes both 4-momentum and 
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charge ill defined. "Sensible mathematics" does not allow us to ignore these infinities "be- 
cause we do not want them" . And even if we ignore the infinities at this stage, they will be 
reappearing again and again in various (practically all) calculations. Quantization based on 
the Fock space simply does not correctly work, although in some sense it must be close to 
a correct theory. 



B. Reducible representation of HOLA 

Now assume that 

Oi(p) = \P)(P\®HP), (230) 

bi(p) = \p)(p\®a(p), (231) 

a 2 (p) = l®a(p), (232) 

b 2 {p) = l®S(p). (233) 

I do not specify at this moment the exact mathematical meanings of \p)(p\ and the hatted 
operators — we will return to it. We obtain 

Pptyufa] = J dpp l t(b l (pya 2 {p) + a 1 (p)6 2 (p) t ) 

dpPfi\p)(p\ ® (a(p) f a(p) + b(p)b(p)^ 

dpPn\p)(p\ ® (a(p) ] a(p) + b(p) ] b(p)j + J dpp tl \p){p\ <g> fc(p,p), 

(234) 

Q[<b.,fa] = Q J dp(b x {py a 2 {p) - Oi(p)6 2 (p) t ) 

= q J dp\p)(p\ <g) (a{p) ] a{p) - b(p)b(p)^ 

= q J dp\p){p\® (a(p) ] d(p) - S(p) f S(p)) - q J dp\p)(p\ <g> t(p,p), (235) 

where t(p,p) = [b(p), b(pY] is the analogue of (I129p . The resulting operators are, of course, 
Hermitian. 

As we can see, the 4-momentum and charge involve number operators of the form (11241) . 
Several approaches to relativistic field theoretic formalism, based on t(p,p) = 1, were de- 



scribed in detail in 



24 



271 ]. Later on in these notes we will reconsider all the steps involved 



in such a choice of i(p,p). 
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C. Sequential approach 



The next application is the following. Assume we have two sequences of operators, a(p, n), 
b(p,n), convergent in some sense to a(p) and b(p). Now take 



\\Jpx 



<h(x) = J dp(a(p,n 2 )e- ipx + b(p,n 2 )^e ipx 
4>i(x) = J dp(b(p, n^)e~ vpx + a(p, n{fe v - 
Then 

f M [0i,02] = J d 3 xT fl0 (x ,x;n 1 ,n 2 ) (236) 

= J dpp^^a(p,n 1 ) j a(p,n 2 ) + b(p,n 1 )b(p,n 2 ) j ^ , (237) 

q[fa,<fa] = J d 3 x J (x ,x;n u n2) (238) 

= qj dp(a(p,n 1 ) f a(p,n 2 ) -b(p,n 1 )b(p,n 2 yy (239) 

are time independent for all n\ and n 2 . We will later consider situations where the Fourier 
form of Dirac delta (1212)1 is rigorously applicable in (I236p and ( 1238jl if n\ and n 2 are finite. 
Integration in (1236 j) and (I238P must be then performed before the limits nj —> oo, n 2 — > 
oo are evaluated. Alternatively, the limits n\ — > oo, n 2 — > oo may be performed under 
momentum integrals in (12371) and (12391) . but not under position integrals in (I236P and f !238[) . 



VIII. SEQUENTIAL APPROACH TO DIRAC DELTAS 

Dirac described his 8(x) as follows: "To get a picture of S(x), take a function of the 
real variable x which vanishes everywhere except inside a small domain, of length e say, 
surrounding the origin x = 0, and which is so large inside this domain that its integral over 
this domain is unity. The exact shape of the function inside this domain does not matter, 
provided there are no unnecessarily wild variations (for example provided the function is 
always of order £ _1 ). Then in the limit e — > this function will go over into S(x)" 28]. 

One possible meaning of Dirac's delta was formalized by Laurent Schwartz [29[], who 
interpreted it as an "evaluation-at-zero map", i.e. a linear functional that, given a function 
/, returns its value at x — 0. In the standard functional-analysis notation we would write 
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(^1/) = /(O) or > more generally, (S x \f) = f(x); in the bra-ket notation of Dirac one would 
expresses the same as (x|/) = f(x). 

An alternative to the Schwartz formalism was developed by the Polish mathematician 
Jan Mikusihski [30J and his coworkers, a program that culminated in the textbook 31] . The 
authors wanted to make precise the intuition of Dirac that 



dx 5(x)f(x) = lim / dx 5 n (x)f(x), (240) 

: J — OO 

and, on this basis, formulate the whole of theory of distributions. In introduction to jsjj 
they wrote: "We shall not avail ourselves to the methods of functional analysis and we shall 
not define distributions as functionals. In applied mathematics distributions are regarded as 
ordinary functions, e.g. the function S(x) of Dirac. Essentially, however, distributions are not 
functions but in an intuitive sense, they may be approximated by functions. Approximation, 
strictly defined, is our starting point for the definition of distributions." 

The idea of "filtering integrals", satisfying lim^oo J"^° dx 8 n (x)f(x) = [/(0_) + /(0+)]/2, 
had been known to Cauchy, Hermite, Poisson, Kirchhoff, and Heaviside many decades before 
Dirac rediscovered the notion for the purposes of quantum mechanics (for a review see the 
fifth chapter of 32| ) . The Mikusihski sequential approach goes deeper, generalizing the 
idea to general distributions. Distributions are in this formalism equivalence classes of 
fundamental sequences of ordinary functions, much the same way as in Cantor's theory real 
numbers are defined as equivalence classes of fundamental sequences of rational numbers. It 
follows that one first has to define fundamental sequences of functions, and then introduce 
an equivalence relation that leads to equivalence classes. 

Let me begin with recalling some basic notions. We say that f n converges uniformly to 
/, if for any e > one can find n such that \ f n (x) — f(x)\ < e for any n > n . The sequence 
fn( x ) = x/n is not uniformly convergent to f(x) = 0, since no matter what e > and n < oo 
we take, we will always find x for which \x/n — 0| > e (just take any x > ne). We say that 
f n converges in the interval A < x < B almost uniformly to /, if it converges to / uniformly 
on each finite closed interval contained in the interval A < x < B (this definition of almost 
uniform convergence is employed in 3jJ). The sequence f n {x) = x/n is almost uniformly 
convergent to f(x) — on the whole real axis since x/n converges uniformly to on each 
interval — oo <a<x<6<oo (now \x\ cannot be arbitrarily large and \x/n — 0| < e, for 
each x in the interval, if n is sufficiently large). 
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In part I of 31] a delta-sequence is defined as any sequence 5 n of continuous functions, 
satisfying 

/oo 
dx8 n (x) = 1, (241) 
-oo 

5 n {x) = (242) 

where F n are twice different iable functions almost uniformly convergent to 

{0 for x < 
(243) 
x for x > 

and the prime denotes the derivative. 

The intuition behind the construction is the following. First, 

{0 for x < 
(244) 
1 for x > 

The derivatives at are defined by right and left limits. 8'(x) is zero everywhere, with the 
exception of x = 0, where 6{x) "jumps from to 1 infinitely fast" — this is, roughly, how 
Dirac imagined his delta "function" . 

A nontrivial and useful example of a delta- sequence is obtained by means of Cauchy's 
principal value, 

• i , 

/ 

-1/e 

If f(x) does not grow too fast with x — > ±oo (i.e. when for large n the function f(x/n) is 
sufficiently slowly changing if compared to 1/x) then 

/oo e ix , r-e l' 1 / e \ e ix 

dxfix/n) — = /(0) lim ( / + / )dx — 
x e-+o+\J_ 1/e J e J x 

POO 

= inf(0) = in / dxf(x)5(x). (246) 



/oo gixn roc ^ix , r—t P^/ e \ ^ix 

dxf(x) = / dxf(x/n) — = lim ( / + / )dxf(x/n) — . (245) 
■oo X J-oo X e ^°+ \J-l/e Je ' X 



In this sense 



— = S n (x). (247) 
mx 

As we can see, delta-sequences can be given also by complex-valued functions. 

In parts II and III of the book the authors further restrict delta- sequences to the narrower 
class of those S n that are smooth, and S n (x) = for \x\ > a n , where a n is a sequence of 
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real numbers convergent to (now 5 n (x) is exactly vanishing outside a given integral - 
the shorter, the greater n). The authors explain in the introduction that there are certain 
differences between different parts of the book since they correspond to different versions of 
the theory, developed by the Mikusihski group over different time periods. In particular, the 
requirement of smoothness and exact vanishing of S n (x) outside of an interval is not needed 
for (1240p to hold, but turns out convenient for some applications. 

I stress the latter property of the Mikusihski-Antosik-Sikorski approach since in what 
follows I want to perform a similar step and restrict admissible delta-sequences to functions 
with some specified properties. 

Now, as we have developed some intuitions for the sequential approach, let us describe in 
more detail the equivalence relation that defines Dirac's delta. We say that the sequence of 
continuous functions f n defined for A < x < B is fundamental if there exists an integer k > 
and an almost uniformly convergent sequence F n of functions, satisfying d k F n (x)/dx k = 
f n {x). k = 2 for S n . The sequences f n and g n are equivalent (we write /„ ~ g n ), if 
d k F n (x)/dx k = f n (x), d k G n (x)/dx k = g n (x), for the same k, and both F n and G n converge 
almost uniformly to the same function. 

One can prove (cf. p. 11 in |3l|) that ~ is an equivalence relation. The equivalence class 
[fn] ='■ f is the distribution associated with equivalent fundamental sequences. Any f n from 
the equivalence class is a representative of the distribution /. 

Let us explicitly show an appropriate construction of 5. Let 5 n be the function of the 
type shown in Fig. 1. It is clear that plots similar to the middle and lower ones would be 
obtained if one replaced this concrete form of 5 n by any function with the same support 
(and integrable to unity), or by a sufficiently narrow Gaussian supported on the whole of R. 
All such delta-sequences thus belong to the same equivalence class 5 = [5 n ] . Delta sequences 
that tend to infinity at x = will be termed A-shaped. 

IX. FURTHER SPLITTING OF EQUIVALENCE CLASSES: DISCONTINUITIES 
AND DELTA-SEQUENCES REGULAR AT ZERO 

In Eq. ( I240p I have purposefully avoided one more natural identification, namely 



•oo 



dx S(x)f(x) = lim 



dxS n (x)f{x) = f{0) = (S\f). 



(248) 



— oo 
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FIG. 1: A-shaped 5 n (x) (upper), 9 n (x) = J_ dx\ 8 n {x\) (middle), and F n {x) = 
fH oc dx2 J® 2 dxi5 n {x\) (lower). 5 n (x) = F^{x). The support of 5 n {x) is here given by the in- 
terval [— ^, with n = 8. F n uniformly converges to F given by (|243|) . 



The reason for this omittance is that (I248[) is true only for functions continuous at x = 0. 
In case of discontinuity the delta-sequence shown at Fig. 1 would imply 

/(0_)+/(0+) 



lim 

n— >oo 



dx5 n (x)f(x) 



(249) 



where /(0±) denotes the left and right limits of f(x) at x — 0. 

But what would have happened had we replaced this concrete delta-sequence by a new 
(fundamental) sequence 5 n (x) = 8 n (x—-)l All the plots from Fig. 1 would be simply shifted 
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to the right by 1/n, so this is again a delta-sequence, but 

/"OO 

lim / dx~8 n (x)f(x) = /(0+). (250) 



Of particular interest are also the following two examples (M-shaped delta sequence, Fig. 2) 

8 n {x) = -Sjx--) +^5 n (x + -) (251) 
2 V nJ 2 V n/ 

and (AA-shaped delta sequence, Fig. 3) 

8 n (x) = U n (x - -) + U n (x+-), (252) 
2 V nJ 2 V n/ 



both yielding 



lim / dx8 n (x)f(x) = Ag=j±Ag±l . (253) 



It is clear that the equivalence relation we have discussed in the previous section can only 
be true if Dirac deltas "act" on continuous /. Schwartz and the other authors go typically 
even further and assume that one deals with functions that are smooth. 

However, in quantum field theory one encounters expressions of the form 

dx8{x)8{x)f{x) (254) 

(or even worse), involving "functions" 8(x)f(x) that are quite far from any form of continu- 
ity. This is one of the sources of the infinities that plague "the so-called quantum electro- 
dynamics" (the phrase of Dirac [l|), especially in loop-amplitude calculations. Intuitively 
one expects that ( 12541) equals 5(0)/(0), whatever this means. It is interesting that such a 
rule can be indeed derived from dimensional regularization techniques for path integrals, 
but 5(0) is then treated as the "infinite quantity J dk/(2Tr)" j^ . 

From a sequential point of view an appropriate calculation could read 

oo 



dx 8(x)8(x)f(x) = lim lim / dx 5 n (x)S m (x)f(x) (255) 

n— >co m— >oo / 

J — oo 

/oo 
dx 8 n (x)8 m (x)f(x) (256) 
-oo 

= lim L(0)/(0) (257) 

m— >oa 

= lim 5 n (0)/(0), (258) 

n— >oo 
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FIG. 2: Analogues of the plots from Fig. 1 but for M-shaped delta-sequence (|25ip . F^(0) = but 
the third derivative is not continuous at x = 0. F n uniformly converges to F given by (|243p . 

where S n (x) and S n (x) are, in principle, different representatives of 5 = [S n ] = [S n ]. It is 
obvious that in order that the expression be well defined one has to require 



lim 4(0)/(0) = lim 5„(0)/(0). 

n— >oo n— >oo 

The examples of delta-sequences we have discussed above would imply 

i imW (o)^x ^:^ , 



lim 5 n (0) 



lim <J n (0)/(0. 

l— >oo 

/(0-) + /(0 4 



x /(0+), 
0x^ + ^. 



(259) 

(260) 
(261) 

(262) 
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FIG. 3: Analogues of the plots from Fig. 1 but for AA-shaped delta-sequence (1252[) . F n uniformly 
converges to F given by (I243p . and all derivatives of F' n vanish at x = 0. 



so that 



dx S(x)S(x)f(x) 



(263) 



becomes ambiguous, unless one restricts delta sequences to a subset of [S n ]. In other words, 
we have to modify the equivalence relation. 



Remark: The proof that 5 2 does not exist, given in |3jj, is based on the assumption that 
the sequential definition of ( I254p should read 



dx 8(x)6(x)f(x) 



lim 

n— >oo 



dx5 n (x) 2 f(x) 



(264) 
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which indeed does not exist for any delta-sequence. It seems that the reason for such a 
choice of definition lies in ambiguities implied by the form of ~. However, we are free to 
change the equivalence relation. A 

One option is the following: Two delta-sequences 8 n and 8 n are equivalent if 8 n ~ 8 n with 
respect to the Mikusihski-Antosik-Sikorski relation discussed above and, in addition, 

lim 8 n (0) = lim 8 n (0) =: 5(0), (265) 

n— ¥00 n— >oo 

/oo /»oo 

dx8 n (x)f(x) = lim / dx8 n (x)f(x) 
~, n— >oo / „ 
-oo J —oo 



= \fiO-) + ^/(0+), (266) 

for any /. Equivalence classes with respect to this new equivalence relation define different 
Dirac deltas. Deltas belonging to the same equivalence class can be uniquely multiplied, 

/oo 
dx8 { n \\x) ...8 { J ) (x)f(x) (267) 
-OO 

= 7(0-) + ^(0 + r 7(0+) (268) 

= ^(0r^(/(0-) + /(0 + )). (269) 

Here 8$ (x) . . . 6$ (x) means that we are free to take any (continuous) representatives 8n](x) 
for any of the delta- sequences. The examples show that there is no difficulty with assuming 
even 8(0) = 0, an option suggested by relativistic invariance, as we shall see in a moment. 



Delta-sequences vanishing at are not a new concept (cf. Chapter V, Eq. (38) in [32j). 
An intriguing example of such a delta-sequence was recently found in the context of time- 
of-arrival operator 
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35| . In spite of this, many authors who generalize the concept of 



Dirac's delta stick to the "obvious" requirement of either divergence or indefiniteness of 8(0) 



cf. jafiHaai). 



X. M-SHAPED AND M-SHAPED DIRAC DELTAS 

In what follows we will mostly work in "momentum space" so that delta-sequences will 
depend on arguments p, k, etc., whereas x will be reserved for their Fourier-transform 
arguments. This is perhaps different from what one is accustomed to, but is more convenient 
from the point of view of applications we have in mind. 
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FIG. 4: The M-shaped function (1270p with a = 1, e = 1/2 (upper), and its Fourier transform 
(lower). 

A. M-shaped delta-sequences 

Let us consider the function shown in the upper part of Fig. 4. It is a particular example, 
for a = 1 and e = 1/2, of 







for 



k<-l 



(f + 2)(f-l) for 





for 



V 



(270) 



(a > 0, e > 0). Its Fourier transform, 
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1 



oo 



8(x,a,e) = — 8(k,a,e)e dk 



oo 



8ea+ (4-ea)cosf . 2 ex 
= 2~2 sm T' ( 2T1 ) 

lim5(x,a,e) = -L (272) 

is the real function shown in the lower part of Fig. 4 . The sequence S n (k,a) = 8(k,a, -), 
with natural n (i.e. e = 1/n), is an example of what I call an M-shaped delta-sequence, which 



is again a particular example of the delta- sequence in the sense of 3Jj. Indeed, 

/oo 
8(k,a,±)dk = 1, (273) 
■oo 

and 

lim f°° f{k)6(k,a,±)dk = /( °' ) t /(0+) - (274) 

71— >00 / ' 
J — OO a 

M-shaped delta-sequences do not have to vanish at 0, 

5(0,a,i) = a, (275) 

for all n, so that 

lim 8(0, a, J) = a. (276) 

n— >oo " 

For each a we deal with a sequence belonging to a different equivalence class. We will now 
show that a = we have encountered before is an important special case. 



B. M-shaped deltas with respect to more general measures 

Let us assume that instead of dp we have to use a measure dfi(p) = p(p)dp, and an 
appropriate delta is needed, 

J dp(p')8(p,p')f(p') = f(p), (277) 

with 8(p,p) = a, say, where a is a constant (that is, not a function of p). The standard 
solution 

8(p,p')= P (p'r 1 5(p-p'), (278) 
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if generalized to M-shaped deltas implies 5(p,p) = p(p)^ 1 5(0) and will not lead to a indepen- 
dent of p (unless 5(0) = a = 0). 

So let us try a different option. Let 5(p,a, ^), 5(0, a, ^) = a, be an arbitrary M-shaped 
delta- sequence. The sequence 

Sn(p, P ') = p(pT 1 Kp-P^p(p)^), (279) 
p(p) = dp(p)/dp, has the required properties 

lim / dp(p')5 n (p,p')f(p') = (280) 
5 n (p,p) = a (281) 
As an exercise, compute in two different ways 
dp(p')5(p,p')5(p + k,p' + k)f(p') 

= lim lim / dp' p(p')p(p')~ 1 5(p - p' , ap(p), ±) p(p' + k)~ 1 5(p - p' , ap(p + k), ^) f(p') 
= lim 5(p - p, ap(p), ^)p(p + k)- 1 f(p) 

n— >oo 

= ap(p)p(p + k)- 1 f(p). (282) 
Reversing the order of limits we would get 
J dp(p')5(p,p')5(p + k,p' + k)f(p') 

= lim lim dp' p(p')p(p'Y l 5{p - p ,ap(p), p(p' + k)~ 1 5(p - p' ,ap(p + k), -M f(p') 

m— >oo n— >co J 

= lim p(p + k)~ 1 5(p-p,ap(p + k),^)f(p) 
= p(p + k)~ 1 ap(p + k)f(p) 

= af(p). (283) 

Consistency of the two calculations requires that 

ap(p)p(p + k)~ l = a (284) 

for all p and k, so that either a ^ and then p(p) = const, or p(p) ^ const and a = 0. 

In non relativistic quantum mechanics p(p) = const so that a can be arbitrary. However, 
relativistic measures involve a nontrivial p and thus we have to assume a = 0. From now 
on, if I write of M-shaped delta- sequences, as opposed to the M-shaped ones, I mean those 
corresponding to a = 0. A particular class of M-shaped delta sequences are the AA-shaped 
sequences whose derivatives of all orders vanish at 0. 
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C. M-shaped versus A-shaped delta-sequences 



Let us discuss in more detail the properties of delta-sequences of the types shown in 
Figs. 1-3. All of them are derived from a single 5 n (k), and can be regarded as particular 
cases of the formula 

1- /. 7\ 1 



Sn(k,j) = -6 n (k - £) + -5 n (- k - £) = 5 n (-k,j), j = 0,l,2,... 
The Fourier transform 



1 f°° 

6 n(x,j) = 7^ dk5 n (k,j)e lkx 



= l:!- f dk5 n (k-i-) e ikx + 1-]- f dk 5 n ( - k - -)e ikx 



CO 

oo 



1 1 r°° • j 

= / dk5Jk)e ikx e ijx/n + 



dk S n (k)e 



-ikXg—ijx/n 



1 



1- 



-5 n (x)e ijx / n + -6 n (x)e- ijx ' n 



(285) 

(286) 
(287) 
(288) 
(289) 



is real. 

It is instructive to pause here for a moment and perform all these calculations for the 
explicit choice of A-shaped delta- sequence from Fig. 1 (now the argument is k and not x), 

for k < - 1 , 

n 2 k + n for — ^ < k < 0, 

-n 2 k + n for < k < \, 

for k > i. 



<J n (fc) 



(290) 



Its Fourier transform 



5 n (x) 



2ir 



1 / X \ ~ Z r, X 

— — sm 2 — 

2%\2nJ 2n 



dk S n (k)e 



ikx 



(291) 



is real, symmetric, bounded, |5„(x)| < l/(27r), and almost uniformly (but not uniformly) 
convergent to l/(2n). 

For all delta- sequences generated from a symmetric S n (k) = S n (k,0), S n (x) = S n (x), we 
can further simplify 



S n (x,j) = 5 n {x)cos—, 

71/ 



lim S n (x,j) = — lim cos — lim / dk5 n (k)e 

n^>oo 2lX n^oo Tl n-^oo 



J — c 



Akx 



2tt' 



(292) 
(293) 
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The latter formulas indicate very clearly that there cannot be much difference between 
A-shaped delta sequences we are accustomed to, and the regular-at-zero M-shaped and 
AA-shaped ones. In fact, they all define the same distribution in the sense of Mikusihski, 
but 5 n (k, 0) and d n (k,j), j > 0, belong to different equivalence classes with respect to our 
modified relation. In particular, products of Dirac deltas S(k,j) = [5 n (k,j)] can be taken 
for all integer j > 0, but 5(k, 0)8(k,j) is ill defined. 

If both /(0_) and /(0+) are finite, we get / dk 8{k, j) N f(k) = for all N = 2, 3, . . . and 
j = 1, 2, . . . . If any of /(0±) is infinite, the expression is still not well defined. This is why 
a modification of Dirac's delta cannot solve all the problems of field quantization. However, 
the situation is much better if we additionally take a nontrivially chosen I(p) occurring in 
HOLA. I will discuss the issue later on in these notes. 

D. 5[f(k)] for M-shaped Dirac deltas 

Let f(k) be a function that vanishes at some ki, i.e. f(ki) = 0, and let 5 n (k) be a delta- 
sequence such that supports of 5 n {k — ki) and 5 n {k — k 2 ) do not overlap, no matter which 
k\ and k 2 , f{k\) = f{k 2 ) = 0, one takes. I will now show how to derive in the context of 
M-shaped deltas the analogue of the standard formula 



Let us first assume that f'{k{) exists. The first trick is to replace, in a neighborhood of a hi, 
the function f(k) by gi{k) = f'(ki)(k — k{) i.e. its tangent line at k\. Then 




(294) 




(295) 



which coincides with (12941) . As an example consider the important case of 



f(ko) 



(296) 



2k , 



(297) 



±Vk 2 + m 2 . 



(298) 
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Let ki = -Vk 2 + m 2 . Then 



gi (ko) = -2Vk 2 + m 2 (k + Vk 2 + m 2 ), (299) 
g 2 (k ) = 2Vk 2 + m 2 (k - Vk 2 + m 2 ). (300) 

For m = the point k = is singular since f'(ki) = f'(k 2 ) = 0. In expressions such as 

<j>(x) = J d 4 k5(k 2 -k 2 )j)(k ,k)e- ikx (301) 

the function (p(k , k) must satisfy 

hm = 0. (302) 

fco->-o± ko 

The fact that 5 is M-shaped is thus irrelevant in this context. 

XI. PLANE WAVES AND M-SHAPED DIRAC DELTAS 

In Dirac's bra-ket notation one encounters formulas such as (k\k') = 2n5(k — k') or 
'h f™oc dk\k)(k\ = 1. In what follows I want to show in what sense they can be interpreted 
in the language of M-shaped delta sequences. 

Let S n (k) = 5 n (k,j) for some positive integer j. We start with convolution of two Mi- 
shap ed delta- sequences, 



= S n *8 m (k) = / dk'5 n (k-k')5 m (k') = 5* mn (k), (303) 

J — oo 

/oo 
dk' S n (k - k')S m (k') =5 n (k), (304) 
-oo 

/oo 
dk'6 n (k-k')5 m (k') = 6 m (k). (305) 
■oo 

The new sequence is again a delta-sequence, 

/oo poo poo 

dkS* nm (k) = dk dk' 5 n (k - k')5 m (k') 
oo J — oo J —oo 

/oo />oo 
dk'5 m (k') / dk5 n (k-k') 
-oo J — oo 

/oo 
dk'5 m (k') = 1, (306) 
-oo 

/oo /*oo />co 

dkf(k)S* nm (k) = lim lim / rfA; / dk' f(k)5 n (k - k')5 m (k') 
■oo J — oo J — oo 

= lim r dkf(k)5 n (k) = fi °- )+ o fi0+) , (307) 
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but 

/oo 
8 n (0 - k')8 m {k')dk! 
-oo 

/oo 
8 n {k')8 m {k')dk' (308) 
-oo 

in general depends on n and m. The other properties are nevertheless analogous to M-shaped 
delta-sequences, 

/oo 
8 n {k')8 m (k')dk' = 8 n (0) = 0, (309) 
■oo 

and 

/oo 
8 n (k')8 n (k')dk' = oo. 
■oo 

Moreover, there exists a sequence a nm , 

lim lim a nm = lim lim a nm = 0, (310) 

m— >oo n— >oo n— >oo m— >oo 

such that 8^ m (k) = for |fc| > a nm . Employing the Fourier transform 

KM = Jk8* nm (k)e ik * = 2n8 n (x)8 m (x) (311) 

we can write 

/oo 
dx8* nm (x)e-^ k - k > 
-oo 

/oo 
dx<^)e^<L(x)e ifc ' z (312) 
-oo 

since 8 n (x) is real. The evaluation map (x\tp) = ip(x), returning the value of Fourier trans- 
form of a given ip, can be used to denote 

(x\k,n) = 2n5 n (x)e ikx = 27r8 n (x)(x\k), (313) 
lim (x\k,n) = e ikx =: (x\k), (314) 

n— >oo 

(k,n\x) = (x\k,n) = lTxl n {x)e- ikx = 2n8 n {x)(k\x) , (315) 
lim(A;,n|x) = e~ ikx =: (k\x). (316) 

rt^oo 

Accordingly 

Km(k-k') = — J dx{k,n\x){x\k',m) (317) 
= ^-(k,n\k',m). (318) 
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FIG. 5: Plots of ()326|) for n = 4 and different values of j: j = (upper) and j = 1 (lower). For 
j = 1 the integral (|327p vanishes due to the oscillating term, whereas the integrand in (|328p is 
non-negative. 

The fact that <5 n (x) is a Fourier transform of a square-integrable S n (k) implies, by 
391 ] . that (a;|fc,n) is square integrable for finite n, 



Plancherel's theorem 

dx \(x\k,n)\ 2 
A similar result holds for 



dx \27rS n (x)e lkx \ 2 = (2vr) 2 / dx \S n (x)\ 2 < oo. (319) 



1 



< 



(27T) 2 
1 

(2^ 



| (k, n\k', m) | 2 

(k , n\k , n) (k' , m\k' , m) < oo. 



Let us finally return to f)312p and take the (well defined) limit, 



(320) 
(321) 



lim 5* nm (k-k') = 2n lim / dx 5 n (x)e ikx 5 m (x)e ik ' x , 

wo ° J-oo 

dx 5 n (x)e- i{h - k ' )x = 5 n (k - k'), (322) 
which is just the inverse Fourier transform of a square integrable function. The next limit 
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would be, however, purely formal 

1 f°° 

lim 5 n {k -k') = — dx e - i( - k - k ' )x = 5{k - k!) = 5*{k- k'), (323) 

n->oo 27T J_ 00 

still showing that the familiar expression for the Fourier form of Dirac's delta is true also 
for deltas constructed by means M-shaped delta sequences vanishing at 0. But, then, what 
about the celebrated divergency 



1 



oc 



5(0) = — = oo? (324) 

Let us do the calculation more precisely, 

/oo 
dx5 n (x) = 5 n (0). (325) 
-oo 

Inserting our explicit example, 

L{x) = k(£) 2sin2 li cos f ' (326) 

we find that the integral vanishes for M-shaped delta sequences (j = 1,2,...) 

— [°° dx( —) ~ 2 sin 2 — cos 3 — = 0, (327) 
2tt \2nJ 2n n ( J 

and for the A-shaped one (j = 0), 

— dx( —) ~ 2 sin 2 — = n, (328) 
2tt 7_ 00 \2nJ 2n ( J 

diverges in the limit n — > oo, as expected (compare Fig. 5). 



XII. M-SHAPED DIRAC DELTAS AND SPECTRAL THEOREM 

Following Dirac, we will denote scalar products of square-integrable functions by the same 
bra-ket symbol as the evaluation map. Let ip(x) = (x\ip) be square-integrable. Then 

/■OO 

(k,n\ip) = / dx (k, n\x)(x\ip) 



■oo 
oo 

dx2n5 n (x)e- ikx tp(x). (329) 



is a kind of windowed Fourier transform 



39[ of ip(x), with the window function 27r5 n (x). 
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Let us again be as explicit as possible and employ the concrete A-shaped sequence ( I290p . 

1 / X \ -2 



(k,n\ip) 



dx 2n- 



— oo 
oo 



2tt \2n 



dx e 



-ikx 



dxe ikx ip n (x,j). 



x 

2n 



sin 2 — cos— e ikx ip(x) 
2n n 

sin — cos — tpix) 
2n n 



(330) 



The Parseval theorem, 
1 

2^ 

with 



dk (ip\k, m)(k, n\ip) 



dx ip m (x,j)i(j n (x,j), 



(331) 



\^m{x,j)i) n (x,j) \ = |</>(x)p 



\2mJ 



sin cos — — 

2m m\2n 

V 



o X 7 X 

sin — cos — 

2n n 



<i 



allows us to use the dominated convergence theorem (Theorem 6 in 40J|) and take the limits, 
in any order, under the integral sign, 



1 f°° 

lim lim — / dk {ip\k , m) {k , n\ 

n->oo m-s>oo 2lT J _ oq 



lim lim 



dx \ib(x) 



\2mJ \2n 



-2 



sm 



— oo 

00 * / T \ -2 / T \ -2 



h/'fa?) lim lim 



oo 



n->oom-too \2lTl/ \2tI 



sm 







X 


JX 




cos — sm 




cos — 


2m 


m 


2n 


n 


X 


• 2 


X 


jx 




cos — sm 




cos — 


2m 


m 


2n" 


n 



dx \ib(x) 



(332) 



The same argument can be applied to the diagonal limit (n = m) 
1 



lim — / dk {ip\k, n){k, n\ip) 

n— >oo 2ll 



oo 
oo 



dx\ijj(x)\ lim 

n— ¥oo 



X \ - 2 

2n 



sm — cos — 
2n n 



dx\ip(x)\ 2 = (?/#> 



— / dk lim (ip\k, n){k, n\ijj). 

2-K J_ OQ n^oo 



Expressing scalar products in terms of norms (the polarization identity 

(f\g) = \(\\f + g\\ 2 -\\f-g\\ 2 + m + ig\\ 2 -i\\f-w 



(333) 
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we can extend the result to pairs of arbitrary square-integrable functions ip and 0, 



lim — 

n— ¥00 2,71 



dk (<f>\k, n)(k, n\ip) 



1 

-y poo 

lim lim — 

m— >oo n— >oo 2,71 



1 

<0|^> 



dfc lim n)(/c, n|^) 

dfc (0|/c, to) (/c, n|-?/>) 
dfc lim lim m)(k, n\ip) 



m—toc n->oo 



(334) 



In this sense 



I = lim — 

n— >oo 2t\ 



dk \k, n)(k, n\ 

> 

i z"°° 

= — / dk lim n)(fc, n| 

i Z"°° 

= lim lim — / dk\k,m)(k,n\ 

m^oo n-s>oo 271 > /_ 00 

i z"°° 

— — I dk lim lim m)(fc, n| 

Z7T J — oo m— »oo n— >oo 

1 Z°° 
= ^ J dk\k)(k\. 



(335) 



Now consider the matrix element 



1 Z°° 

((f>\A\ij;) = lim lim — / dk A(k){(f)\k, m)(k, n\ip) 

m->oo n->-oc 2% > /_ 00 

OO POO 

dk A(k) / c&c (x\k, to) I dy (k : n\y)(y\ip) 



= lim lim — 

m^oo n— >oo 27T 



lim lim — / dk A(k) / dx 4>(x)27i5 m (x)e ikx / dy2n5 n (y)^(y) 

rwoo n->oo 27T J _ oc 7-oc J-oo 



-iky 



OO 
OO 



j /»00 /"OO /"OO 

lim — / dkA(k) / 27r5 m (x)0(x)e- ifea: / dy 27iS n (y)^(y)e- iky . 

>n^oo27lj_ 00 7-cc 7-oc 

With our definitions of Fourier transforms, 



= lim lim — 



/(*) 
/(*) 

dx f(x)g(x)e~ lkx 



1 

27 



2vr 7_ 



rfA;/(A;)e ifcx , 
dxf(x)e~ ikx , 

dk'f(k-k')g(k') : 



(336) 
(337) 
(338) 



we obtain 



-I />oo />oo roo 

= lim lim — / dkA(k) dk' S m (k - k')^{F) dk" S n (k - k")^{k") . 

m— >oo n— »oo 27T / ~, / ~, / 

" J —oo J — oo «/ — oo 



<5 m *0(fc) 



<5 ra *i/>(fc) 
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As a by-product of the above calculation we have shown that (k,n\if)) = S n * ip{k). 
Of great importance is the following 



Lemma: (a variant of Theorem II. 3. 1.1 in 
converges to ip(k) almost uniformly. 



311 ] ) If ip{k) is continuous then (k,n\ip) = S n *ip(k) 



Proof: Assume for simplicity that the delta-sequence is non-negative, S n (k) > 0, and van- 
ishes for \k\ > a n , linin-s.oo a n = 0, like in all our previous examples. Assume k belongs to a 
finite interval [a, b]. Then 



\6 n * if)(k) - ip(k)\ = 


poo 

/ dk' 5 n (k 

J — oo 


- k')^(k') 




POO 

1 dk' 5 n (k — k') 

—oo 




POO 

/ dk' 8 n (k 

J — oo 




- i>{k) 


) 


<-) 


POO 

dk' 5 n (k - 

—oo 


■ k')\ip(k') - 






= / 


fOt n 

dk' 8 n {k') 


\ip(k - k') 


-*i>(k)\ 





By definition of continuity of ip(k) at k, for any e > there exists no such that \ip(k — k') — 
i^{k)\ < e for all n > n and — a n <k'< a n . So 

\5 n *ip(k)-ip(k)\ < / " dk' 8 n (k')\ip(k - k') - ip(k)\ 



dk' 8 n {k') = e 



< e 



for all n> tiq and a < k < b.M 

We can take the limits under the integral sign in 

1 



lim lim 

m— too n—too 2,1T 



dk A(k)5 m * <p{k)5 n * 4>(k) 



in several cases. The simplest is the one of A(k)5 m * <p{k)8 n * 4>{k) converging uniformly to 



A(k)5 m *(J)(k)i/j(k) and A(k)(f)(k)5 n *i/j(k) . For, example let A(k) be continuous for a < k < b, 
and A(k) = for k > b, k < a. Almost uniform convergence implies uniform convergence 
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on any closed interval, and 



1 r b 



(<f>\A\ij;) = lim lim — / dk A(k)5 m * <p{k)5 n * ip(k) 

m— >oo 71— >oo 2ll J „ 
1 f b 

= — / dk A{k) lim S m * </>(k) lim 8 n *ip{k) 

2tT l„ m— >oo 71— >oo 



1 ' b 



— / dkA(k)<f>(kWk) 
2ii 



1 

2^ 



b 

dkA(k){(j)\k)(k\^). 



a 



A generalization going beyond finite intervals a < k < b can be based on the observation 
that dominated convergence theorem guarantees 

-I /*00 1 /"OO 

lim lim — / dkA(k)5 m * fi(k)5 n * ip(k) = — dk A{k)J{k)i){k) (339) 

m->oo n->oo 2lT J _ QO 27T J _ !XJ 

if the right-hand-side of (13391) is finite and 



\A(k)5 m * (j)(k)8 n * ip(k)\ < consti x \A(k)5 m * ^(k)i/j(k)\ < const 2 x \A(k)</>(k)ip(k)\. 

Let us now turn to the issue of eigenvalues of A. Let <f>(k) be continuous. The eigenvalue 
problem for A should be understood in this formalism in the following sense 

(cj>\A\p) = lim(0L4b,O 

1 f°° 

= lim lim lim — / dk A( y k)((f>\k,m)(k,n\p,l) 

l— >oo m— >oo n— »oo 2/T I ^, 

/oo 
dk A(k)((f)\k,m)S n i(k -p) (340) 
-oo 

= lim \(A{p-) + A(p + )) (<f>\p, m) 

771— >O0 / \ / 



= l -(A{ V _) + A{ V+ )){<P\ V ). (341) 

Let us recall that (<fr\k,m) converges almost uniformly to a continuous ((fi\k) = <f)(k), while 
the support of 5 n i(k — p) is compact. Therefore, the sequence A(k)((j)\k, m)5 n i{k, — p) con- 
verges uniformly to A(k)((p\k)5 n i(k — p) if A(k) is bounded on any compact subset of R. 
In consequence, one can alternatively compute the eigenvalue problem as follows [compare 

in; 



oo 



(4>\A\p) = lim lim / dkA(k)(<p\k)5 n i(k- V ) (342) 

l— >oo n— >oo J _ 

= l -{A{p^) + A(p + )){<P\p). (343) 
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If continuous 4>(k) is itself of compact support then (<f>\k,m) converges to (4>\k) uniformly. 

Although f)34ip shows that the sequential approach automatically defines eigenvalues 
also at points of discontinuity, in what follows I simplify discussion and assume that A(k) 
is continuous. In such a case the formalism based on M-shaped deltas does not seem to 
essentially differ from other mathematical approaches to the Dirac bra-ket formalism for 
continuous spectra (cf. 4114431] ). Anyway, what we modify is the value of (k\k') on the 
diagonal k — k' . Looking at the rigged Hilbert space approach, say, one does not find any 
point where the formalism employs a concrete value (infinite or not) of (k\k). 

The main difference with respect to more standard formalisms thus lies in definition of 
diagonal elements of operators whose spectrum is continuous. Indeed, 

(p|v4|p) = lim lim {p, r\A\p, s) 

r— >oo s— >oo 

1 f°° 

= lim lim lim lim — / dkA(k)(p,r\k,m)(k,n\p,s) 

r— >oo s— >oo m— ¥00 n— ¥00 2,71 I 

J — oo 



-oo 



/oo 
dk A(k)S rm (p — k)S ns (k — p) 
-oo 

lim lim 2vr A(p)5 ns (0) = 0. (344) 



s— >oo n— >oo 



In particular, (k\k) = 0. The latter, of course, does not mean that \k) = lim^oo \k,n) is 
vanishing. Simply, as has been stressed here many times, neither (k\k') nor (k\ip) should be 
regarded as scalar products. 

It follows that \k)(k\ is not a projector, 

\k)(k\k)(k\ = 5{0)\k)(k\ = 0. (345) 

However, let Xx(k) be the characteristic function of X C M, i.e. Xx(k) = 1 if k 6 X, and 
Xx(k) = if k X. The operators 

E{X) = -L f°° dk X x(k)\k)(k\, (346) 

^ J-oo 



E([a,b]) = ±- I dk\k)(k\. (317) 

Z7T 



a 
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are projectors (0 is the empty set): 



E(X)E{Y) = E(XnY), (348) 

£(0) = 0, (349) 

E{R) = I, (350) 

i r°° i r b 

E (fab])7T- / dkA(k)\k)(k\ = — / dkA(k) \k)(k\. (351) 

27T 27T ./„ 



|l2| 

(^ifl^WW) = lim ... lim -f- / dpip(p)(p, ni\k,n 2 )(k,n 3 \-^- I dp ip{p')\p ,n 4 ) 



The operator is not a projector but a POVM [17], since 



"J— ^ << 1— X) 27T ./_oo 27T 



—00 
00 



1 r° ,1 



lim ... lim — / dp— / dp' i/j(p)ip(p')(p, rii\k, 712) (k,n 3 \p' ,114) 



/oo /*oo 
dj9 / dp' 1p(p)lp(p')5 nin2 (p - k)5 n3n4 (k - p') 
■00 J-00 

= (^|A;>(A#> = |^)| 2 > 0, (352) 
and I = fdk/(2ir)\k)(k\. 

Remark: The above calculations show what should be meant by products of operators. For 
example, let A, B, C commute. Then 

ABC = I ^A(k)\k)(k\ ^-B(k')\k')(k'\ —C(k")\k")(k"\ 

OO t/ — OO «/ — OO ' 



/°° rlk f°° rlk' 

—A{k)\k,n 1 )(k,n 2 \ / —B{k')\k',n' 1 )(k',n! 2 \ 
-00 27r ' J -00 27r 

^—C(k") \k", n'[) (k", n'zl (353) 

.~ 27T 



= lim ... lim / — — A(k)B(k')C(k'')\k,n 1 ){k,n 2 \k',n' 1 ){k',n' 2 \k'',n'l){k'',n'^\ 
^A(*)B(*)C(*)|*>(fc| 

-00 2?r 

independently of the order of limits, the latter property being equivalent to associativity of 
the product. A 

As we can see, one can work with M-shaped deltas with no difficulty Restricting the 
class of Mikusihski deltas to a subset of [5 n ] we obtain the usual Dirac deltas but equipped 
with new and useful properties. The fact that we work with kets of "zero length" is no more 
paradoxical than what we face when we deal with those of "infinite length" . 



61 



Paraphrasing physics textbook presentations of S(x) we could say that M-shaped delta 
is a "function which is zero everywhere, but whose integral over any interval containing 
equals unity". Such a "definition" is neither worse nor better than the usual one with 
infinity at — the difference between them is on the set of zero measure! 

Let me end this section with remarks on relativistic M-shaped Dirac deltas. In what 
follows two types of expressions involving Dirac deltas will be encountered. First, 

5< 3 \k) = 5*(h)5*(k 2 )5*(k 3 ) = [Cn 1 (^l)][C 2 n 2 (^)][C 3 n3(^)], (354) 

1 = f d 3 kS* {3 \k), (355) 



is M-shaped in ID and 3D: 8*(0) = 0, <5*( 3 )(0) = 0. Sometimes it may be useful to employ 
the integral formula (I212p for 5^ (p) (under appropriate integrals) but then 5^ (0) is mean- 
ingless. From now on I will stick to the convention that p denotes the spacelike part of the 
4- vector p = (a/j? 2 + m 2 , p) = (po,p), whereas k stands for an analogous part of the null 
4-vector k = (\k\, k) = (fc , k). Relativistically covariant kets will be normalized as follows, 

(p\p') = (2k) 3 2^ p 2 + m 2 5< 3 \p-p') =: S m (P,P'), (356) 
(k\k') = (27r) 3 2|fc| 5* {3) {k - k') =: 5 {k, k'), (357) 
5 m (P,P) = 5 (k,k) = 0. (358) 



The resolutions of unity are 

d 3 k 

f d 3 p f 

/ /o r~, — =?1p)(p1 = / d P\p)(p\ = I 

Jr3 (2tt) 3 2a/p 2 + m 2 Jrz 



r d 



k)(k\ = [ dk\k)(k\ = 1, (359) 

(360) 



Here dp and dk denote relativistic measures on the hyperboloid p 2 = m 2 and the light cone 
k 2 = 0, respectively. 



XIII. ALGEBRA OF FREE FIELD OPERATORS REVISITED 

The ultimate goal of these lecture notes is to discuss a theory of interacting quantum 
relativistic scalar fields: A charged massive boson ip(x) and a neutral massless boson (f>(x). 
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For free fields one finds 

i){x) = j ' dp(a(p)e- ipx + b( P y e ipx y (361) 
4>{x) = j ' dk[c{k)e- lkx + c(fc) t e fc ). (362) 
The amplitude particle and antiparticle operators are assumed to satisfy bosonic HOLAs: 

[a(p),a( P y] = 8 m (p,p')I m (p), (363) 

[a(p),n + (p')} = 8 m (p,p')a(p), (364) 

[a(p)\n + (p')} = -5 m (p,p')a(p)\ (365) 

[b(p),b(p'y\ = 5 m (p,p')I m (p), (366) 

[b(p),n-tf)] = S m (p,p')b(p), (367) 

[b{p)\n_{p>)\ = -5 m (p,p')b( P y. (368) 

The massless neutral field satisfies a similar algebra but with the light cone deltas, 

[c(fe),c(fe / ) t ] = 5 (k,k')I (k), (369) 
[c(k),n (k')] = 5 (k,k')c(k), (370) 
[c(k)\n (k')} = -5 (k,k')c(k)l (371) 

The remaining commutators vanish. 

Remark: Let us immediately note the peculiarity of M-shaped deltas: 

[a(p),a( P y\ = [b(p),b(p)t] = [c(fc),c(fc)t] = 0. (372) 

This type of commutation relations is, in fact, implicitly employed in standard quantum 
field theory whenever one applies the normal ordering operation (: q : instead of q, etc.). I 
think this is yet another argument for quantization in terms of M-shaped deltas. A 

Remark: Let us divide R 3 into disjoined sets X n) IJ^Lo^™ = ^ 3 > anc ^ ^ Xn(p) De their 
characteristic functions satisfying Y^=oXn{p) = 1 and Xn(p)Xj(p) — Xn(p)8nj- Here 5 n j is 
the Kronecker delta. Defining 

an = J dpxn(p)a(p) (373) 



63 



we find that 



dpXn{p)a{p), J dp Xn{p')a(p') j 
dp / dp Xn(p)Xn(p)Hp),a(p') j ] 



dp J dp' Xn(p)Xnip')S m (p,p')I m (p) 

dpXn{p) 2 I m {p) 

dpXn(p)Irn(p) = I m ,n ^ 0. 



(374) 



The example shows that commutability f !3T2j) typical of M-shaped deltas does not imply 
that a discrete version of the representation should be automatically taken in the form 
[a n , ot] = 0. The discrete counterpart of HOLA will thus read 



[a n , a) 



(375) 



with appropriately chosen I n . ▲ 

The above fields will be used by me in two nontrivial toy models. The first model involves 
interaction term (in interaction picture) 



Hi(t) 



d 3 xj (t, x)(f)(t,x) 



where 



(376) 

iq I d 3 x (if)(t, x)tp(t, x) - ip(t, x)ip(t, x)^j <j>(t, x) (377) 
d z xp{t } x)(t){t,x), (378) 



(379) 



jn(x) = iqip(x)d fl ifj(x) - iqd^x^x) 



is the charge-density current of the massive field, and p(t, x) is the charge- density operator, 



/ d 3 xp(t,x) = / d 3 xp(0, x) = q[ip\ ip). 

JR3 J^3 



(380) 



The fields that define the Hamiltonian are free, a fact typical of the interaction picture. 
Note that I do not assume normal ordering of p(t, x). 

The Hamiltonian has the essential properties of quantum-electrodynamical interaction, 
but is not relativistically covariant. For this reason I will also consider a more complicated 
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interaction, with the same charge current, 

i?i(t) = J^(x)j^t,x)(t>(t,x), (381) 

where E is a hypersurface in Minkowski space and da^(x) a measure multiplied by a vector 
field of vectors normal to E. Analysis of (at least) loop diagrams for the latter interaction 
is the principal goal of these lecture notes. 



A. Reducible N > 1 representation of HOLA — an indefinite-frequency harmonic 
oscillator representation of ip(x) 

Let [a(p), a(p'Y] = i + (p,p'), [6(p), b(p'Y] = i-(p,p'), where i± are any functions (or 
distributions) satisfying t±(p,p) = 1. We assume that all commutators between the as 
and the 6s vanish. Then N = 1 reducible representation is constructed analogously to the 
indefinite- frequency harmonic oscillator discussed earlier in these notes. The particle part, 

a(p,l) = \p)(p\®a(p), (382) 

a(p,iy = \p)(p\®a(p)\ (383) 

n+(p,l) = |p)(p|®a(p) t a(p), (384) 

I m (p,l) = |p)(p|®l, (385) 

dpl m (p,l) = I m (i), (resolution of identity) (386) 

and the antiparticle part, 

6(P,1) = |p)(p|®S(p), (387) 

6(p,l)t = |p)(p|(g)6(p)t, (388) 

n_(p,l) = |p)(p|<8)6(p)t6(p), (389) 

/ m (p,l) = |p>(p| <8) 1, (390) 

can be explicitly constructed in various ways. The simplest case occurs if a(p), b(p) are 
p-independent, i.e. 

a(p,l) = \p){p\®(a®l), (391) 

a(p,l) f = |p)(p|®(a t (g)l), (392) 

I m (p,l) = |P>{P|®(1® 1), (393) 

n+(p, 1) = |p>(p| (8) (a f a (8)1), (394) 
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6(p,l) = |p)(p|® (l<8>a), (395) 

fe(p,l)t = |p)(p|<8)(l(g)a t ), (396) 

7 m (p,l) = |p)(p|®(l®l), (397) 

ra_(p,l) = |p)(p| <g> (1 (gia^), (398) 

with [a, a'] = 1. At the other extreme is L±(p,p') that equals an M-shaped Dirac delta 
normalized to 1 at 0. Both cases are worthy of consideration. 
For iV > 1 

a(p, N) = J= (a(p, 1) <8> 7 m (l) <g> • • • <g> 7 m (l) + • • • + J m (l) <g> • • • <8> 7 m (l) <g> a(p, 1)) , 

(399) 

a(p, A) f = -j= (a(p, l) f (8> 7 m (l) <8> • • • <8> 7 m (l) + • • • + 7 m (l) <8> • • • <8> 7 m (l) <8> a(p, l) f ) , 

(400) 

n+(p, AO = n+(p, 1) <g> 7 m (l) <g> • • • <8> 7 ro (l) + • • • + 7 m (l) <g> • • • <8> 7 m (l) <g> n+(p, 1), 

(401) 

7 m (p, N) — — (l m {p, 1) ® 7 m (l) <8) ... (8) 7 m (l) + • • • + 7 m (l) <8> • • • <8> 7 m (l) <8> 7 m (p, 1)), 

(402) 

I m (N) = JdpI m (p,N), (resolution of identity) (403) 



b(p, N) = — (b(p, 1) <8> 7 m (l) <8> • • • (8) 7 m (l) + • • • + 7 m (l) (8) ... ® 7 m (l) <8> 6(p, 1)) , 

(40< 

b(p, A)t = _L (fc(p, i)t ® / m (i) ® . . . ® / m (i) + . . . + / m (i) (g, . . . (g, / m (i) ® 6( p> i)t j , 

(40! 

n_(p,A) = n_(p,l)(8)7 m (l)(8)--.(8)7 m (l) + --- + 7 m (l)(8)-..(8)7 m (l)(8)n_(p,l). 



(406) 



Note that 



a(p, l) 2 = \p) (p\p) (p| ® a(p) 2 = (407) 
o 

and similarly with products of all the other elements of HOLA provided they are taken at 
p' = p. In particular, 

n+(p,l) = (|p)(p| ®a(p) t )(l®a(p)) = a(p, l) f (l ® a(p)) ^ a(p, l) f a(p, 1) = 0. 
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This does not mean that the algebra is trivial or Abelian. Indeed, the 4-momentum (|234[) 

P M (1) = J dp Pfl \p){p\ ® (a(p) t a(p) +b{p)b{ P y 

= J dpp tl \p){p\® (a{p) ] a{p) + 6(p) f 6(p)) + ^ dpp M |p)(p| ® i-(p,p), 

= J dpp^(n+{p,l)+n-{p,l)) + f dppplmip, 1), (408) 
P M (iV) = P M (1) ® 7 m (l) ® ■ • • <g> J m (l) + ■ • • + J m (l) ® ■ • • ® 7 m (l) <g> P M (1), (409) 

possesses the required properties of generator of 4-translations: For any N > 1 



e iP{N)x a(p,N)e 


-iP{N)x 


= a(p,N)e- ix ?, 


(410) 


e iP ( N H(p, N)e 


-iP(N)x 


= b(p,N)e- ixp , 


(411) 


e lP{N)x a{p,N)^e 


-iP(N)x 


= a(p,Ny e ixp , 


(412) 


e iP{N)x b(p,N)^e 


-iP(N)x 


= b{p,Ny e ixp , 


(413) 


e lP{N))x ^(0,N)e 


-iP(N)x 


= iI>(x,N), 


(414) 


e- iP{N)y il>(x,N)< 


JP(N)y 


= ip(x-y,N). 


(415) 



Here ip(x,N) denotes ip{x) but taken in the reducible N >1 representation of HOLA. The 
vacuum term 

vacuum N J dpp^I m (p,N) (416) 

is well defined and commutes with all elements of HOLA. In analogy to the standard Fock- 
space formalism we can work with 

:P M (iV) : = P^N) - P M (iV) vacuum = J d PPfl (n + (p,N) +n_(p,JV)) (417) 

but now, as opposed to the standard approach, the subtraction is well defined. 

Remark: We know that one is not allowed to divide by 0. Why? Because then arithmetics 
would become ambiguous. Think of 1 ■ = 2 • 0. If division by would be acceptable then 
we would conclude that 1 = 2. This is the actual reason why 1/0 is not allowed. The 
same situation occurs with 1 + oo = 2 + oo. If we could subtract infinities then 1 = 2. So 
subtraction of infinities is as forbidden in arithmetics as division by 0.A 
Remark: Due to the property of M-shaped deltas, 



J dp Vil (a(p, l)t a (p, 1) + b(p, l)6(p, l) f ) = 
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so this would not be a correct definition of free-field 4-momentum. For higher A 

J dpp, (a(p, Ny a (p, A) + b(p, N)b(p, A)t) 

= J d PP , (a(p,Ny a (p,N) + b(p,Nyb(p,N^) ^0 

but for reasons mentioned earlier I do not think this is the correct free-field Hamiltonian. 
Although, who knows... A 

B. Reducible N > 1 representation of HOLA — an indefinite-frequency harmonic 
oscillator representation of <p(x) 

The constant in the Noether current is C = 1/2. Take c(k) = a. For A = 1 



c(k,l) = 


\k)(k\$ 


3 a, 


(418) 


c(fc,l) t = 


\k)(k\S 




(419) 


7 (fc,l) = 


\k)(k\S 


n, 


(420) 


n (k,l) = 


\k)(k\S 


§ a) a. 


(421) 



For A > 1, 

c(k,N) = -U(c(fc, l)®I (l)®---®Io(l) + --- + I (l)®---®I (l)®c(k, 1)Y 
v A ^ / 

(422) 

c(fc, A) f = -U (c(p, l) f <g> J (l) ® • • • ® 7 (1) + • • • + 7 (1) ® ■ ■ ■ ® / (1) ® c(fc, l) 1 ) , 
v A ^ / 

(423) 

n (fc, A) = no(fe, 1) <E> io(l) ® • • • ® I (l) + ■■■ + 7 (1) ® • • • ® 7 (1) ® n (fe, 1), 

(424) 

7 (fc, A) = 1 (/„(*!, 1) ® 7 (1) ® • • • ® 7 (1) + • • • + 7 (1) (8) ... ® 7 (1) ® 7 (fc, 1)) , 

(425) 

7 (A) = JdkI (k,N), (resolution of identity) (426) 
The indefinite-frequency-type 4-momentum 

P M (1) = J dkk^\k)(k\®^d+^ J dk k^ \k)(k\ ® 1 

= J dkk^n (k,l) + ^ J dkk^Ioik,!) (427) 
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acts generator of 4-translations 

e iP{1)x c(k, l)e- iP{l)x = c(k,l)e- ix \ (428) 

iP(l)x (u i\j -iP(l)x _ „(u ~\\} p ixk 



e w c 



(k, l)t e -^(i> = c {k,iy e ixk , (429) 



e tP(1) >(0, l)e" iP(1)a; = 0(x,l), (430) 
e- lP(1)s/ 0(a;, l)e iP{1)?/ = (431) 



C. Physical interpretation of field operators 

Returning to canonical position operator (I149P 

Q(t, N) = V x j^(a w {N)e-^ + MA^) , (432) 
z — ' V moo V / 

w 

and comparing it with 4>(x, N) evaluated at the origin x = 0, 

M°.*0 = / ,o PES ( c ( fc ' ^)e-^ (fc)t + c(fc, iV)t e -( fc ) f ) , (433) 
J R 3 {2iry2uj{k) V / 

one understands that for N = 1 the field at the origin is a canonical position of some 

indefinite-frequency operator. The field taken at an arbitrary point x is obtained from the 

one at the origin by translation. An elementary quantum field <p(x, 1), corresponding to iV = 

1, consists of a single oscillator that exists "everywhere" in space, i.e. is in superposition of 

different localizations. Representations characterized by iV > 1 have physical interpretation 

of a gas consisting of iV indefinite-frequency noninteracting bosonic oscillators. 

D. (N, iV)-oscillator representation of HOLA: N massive charged and N' massless 
neutral oscillators 

Let % m {N) and %q{N) be the Hilbert spaces of the above two types of N > 1 scalar field 
representations. The full Hilbert space appropriate for the model of interacting fields is 

n(N, N') = n m (N) ® Ho(N') (434) 

with two in principle independent parameters A and N'. The representations are constructed 
by the embeddings 

a(p, N, N') = a(p, N) <g> I (N'), (435) 
c(k,N,N') = I m (N)®c{k,N'), (436) 

69 



and so on. 



XIV. STATES AND THEIR BOUNDARY CONDITIONS FOR (N, N')- 
REPRESENTATION 

We have not yet needed explicit forms of analogues of vacuum, n-particle or coherent 
states of quantum fields. Their construction is analogous to what we have done with states 
of N indefinite-frequency oscillators. 

Vacuum appropriate for representations fl39ip - (l398l) . fl418p -( l421|) . is constructed as fol- 
lows: 

\O m ,l) = J dpO m {p)\p,0,0), (437) 

|O ,l) = J dkO (k)\k,0), (438) 
\O t N,N') = |O m ,l) 07V ® \O A)® N '. (439) 

Vacuum states are normalizable 

J dp\O m (p)\ 2 = J dk\O (k)\ 2 = 1 (440) 

Square integrability implies, roughly speaking, vanishing of O m (p) and Oo{k) in the lim- 
its \p\ — > oo, \k\ — > oo. This is one of those features of field quantization in (N,N')~ 
representations that will in the end make mathematics of field quantization "sensible" (in 
Dirac's sense). For massless fields, such as Oo(k), a boundary condition has to be imposed 
also at k — (ko, k) = (\k\, k) = 0. For massive fields, with p 2 = m? > 0, the point p = 
does not belong to the mass-m hyperboloid: (i^/m 2 + p 2 , p) — > (±m, 0) with |p| — > 0. 

For m = the point k = is Lorentz invariant. It corresponds to k = and k' = Lk = 
L0 = for any L so that k' = as well. The boundary condition Oo(0) = is also Lorentz 
invariant (in fact, it is invariant under the whole Poincare group). Well-definiteness of field 
theory that involves a massless field will require a yet stronger boundary condition, namely 



lim = 0, (441) 



I 



12 



(442) 



for all natural n. The fact that the tip k = of the light cone is Lorentz invariant has 
deep consequences for the structure of unitary representations of the Poincare group and 
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geometry of the light cone. Hilbert spaces of massless k = and k 7^ representations 



are completely different (the representations are induced from different little groups 44j ) 
and there is no Lorentz transformation that could map k = into k! ^ 0, even if k' 2 = 0. 
The light cone splits in a Lorentz invariant way into the tip k Q = and the physical cone 
(ko 0). The requirement that the probability density |Oo(fc)| 2 tends to as k — > is thus 
mathematically very natural. From a physical point of view it eliminates massless fields of 
infinite wavelength and zero frequency. 

The rest of the construction is completely analogous to what I did for iV-oscillator states. 
(However, we have to remember that (p\p) = (k\k) = 0.) In particular, the displacement 
operator is defined as 

V(a,(3,~f,N,N') = exp J dp(a(p)a(p, N, N') j + 0(p)b(p, N, N') j - H.c. j 

xexp J dk(j{k)c{k,N,N'y - H.c.) (443) 

= exp J dp(a(p)a(p,Ny + p(p)b(p,Ny - H.c.) 

g> exp J dk(^r{k)c{k,N'y - H.c.) (444) 

Corresponding statistics of excitations is Renyi-deformed Poissonian. 

XV. RELATIVISTIC COVARIANCE 

Let x be a point in Minkowski space, i.e. 

(x°, x 1 , x 2 , x 3 ) = (ct, x, y, z) = (x , -x x , -x 2 , -x 3 ). (445) 
In what follows we work in units with c = 1, so we set x® = t. The Poincare transformation 

x' = Lx + y, (446) 
where L is a Lorentz transformation and y a 4-vector, is equivalent to 



x 



x' 



L a p x p + y a = x^L- 1 ^ + y a , (447) 
Ljxp + y a = x p L~ l \ + y a . (448) 
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Poincare transformations are represented unitarily by operators U (L, y, N, N') constructed 
below. To begin with, the unitary representation of 4-translations is generated by 

P^(N,N') = J dpp lx (n + (p,N)+n_(p,N)' S )®I (N') + I in (N)® J dkk^n {k,N') 

+N j dpp, I m (p, N) ® Io(N') + I m (N) ® ^ Jdkk, I (k, N') (449) 

A. Action of the Poincare group on field operators 

A pure 4-translation x' — x + y is represented by 

17(1, y, N, N') = e *" W"') = e n"W). ( 450 ) 

It multiplies annihilation operators by phase factors 

U(l, y, N, N'y a (p, N, N')U(1, y, N, N') = a(p, N, N')e™, (451) 
U(l,y,N,N'yb(p,N,N')U(l,y,N,N') = b(p, N, N')e ipy , (452) 
U(l,y, N, Ny C (k, N, N')U(l,y, N, N') = c(k, N, N')e iky . (453) 

Recall that py = a/p 2 + m 2 y°—p-y, ky = \k\y°—k-y. Now let L be a Lorentz transformation 
and Lp be the spacelike part of Lp. Unitary spin-0 representation of the Lorentz group is 
given by 

U m (L,0,l) = J dp\p)(L- x p\®(l®l), (454) 

C/o(L,0,l) = J dk\k){L- x k\®l, (455) 
U(L, 0, N, N') = U m (L, 0, 1)® N ® U (L, 0, 1)®"'. (456) 

U(L, 0, N, N') changes momenta of the amplitude operators (the Doppler effect), 

U (7, 0, N, N')^a(p, N, N')U(L, 0, N, N') = a(L~ 1 p, N, N'), (457) 
U (L, 0, N, N')^b(p, N, N')U(L, 0, N, N') = b(L _1 p, N, TV'), (458) 
U(L,0,N,Ny C (k,N,N')U(L,0,N,N') = c(X _1 fc, N, N'). (459) 

Four-translations and Lorentz transformations are combined into full Poincare transforma- 
tions (L,y) = (l,y)(L,0): 

(L,y)x = (l,y)(L,0)x=(l,y)(Lx + 0) = l(Lx + 0) + y = Lx + y. (460) 
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Therefore 

U(L, y, N, N') = U(l, y, N, N')U(L, 0, N, N') (461) 

is a unitary representation of a general Poincare transformation. As an example of explicit 
Poincare transformation consider the massive charged field 

i[>(x, N, N') = J dp (a(p, N, N')e~ ipx + b(p, N, N')^e ipx ^j . (462) 

Then 

U(L, y, N , N')^(x, N , N')U(L, y, N, TV') 

= J dp(u(L, 0, N, N'yU(l, y, N, N'^aip, N, N')U(l,y, N, N')U(L, 0, N, N')e~ ipx 

+ U(L, 0, N, N'yU(l, y, N, N')%p, N, N')^U{1, y, N, N')U(L, 0, N, N')e ipx ) 

dp(u(L, 0, N, N'y a (p, N, N')U(L, 0, N, N')e~ ip( - x - y) 

+ U(L, 0, N, N')*b(p, N, N'yU(L, 0, N, N')e ip{x ' y) ^ 

= J dp(a{lT x p,N, N')e~ ip(x ~ y) +b(L~ 1 p,N, N')^e ip{x ~ y) ^ 

= J dpfaiL^p, N, N ') e -^L-M--y) + b^L^p, N, N')* e iLL ~ 1 * x - v)> ) 

= J dp(a(p, N, N')e- iLp{x - y) + b(p, N, TV') 1 e iLp{x - y) ^ 

= J dp(a(p, N, N')e- ipL ~ 1{ - x -ri + b(p, N, jv')t e * i_1 (*-tO ) 

= iIj(L~ 1 (x -y),N,N'). (463) 

I have used here the Lorentz invariance : d(L~ l p) = dp and Lp- Lx = p-x = implying 

Lp-(x-y) = Lp- LL~ l (x -y)=p-L~ 1 (x- y). 

L~ x {x — y) is the inverse of Lx + y. Indeed, 

L- 1 ((Lx + y)-y) = x, (464) 
L(L-\x-y))+y = x. (465) 

The same transformation rule holds for the massless neutral field 

U(L, y, N, N'y<f>(x, N, N')U(L, y, N, N') = ^(L'^x -y),N, N') . (466) 
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Number operators transform as translation invariant scalar fields in momentum space 

U(L,y,N,N^n ± (p,N,N')U(L,y,N,N') = n±(L~ 1 p, N, N'), (467) 
U(L,y,N,Nyn (k,N,N')U(L,y,N,N') = n^L^k, N, N'). (468) 

In consequence, the 4-momentum is a 4-vector: 

j dp Pll n±{L- x p,N,N') = J dp (LL~ 1 p) IJ ,n±(L~ 1 p, N, N') 

= j dp(Lp),n ± (p,N,N') 

= L/ J dpp u n±(p,N,N') 
J dkk tl n {L- 1 k,N,N') = J dk (LL- 1 k) f ,n (L- 1 k, N, N') 

= V J dkk u n (k,N,N'), 

so that 

U(L,y,N,N'yP,(N,N')U(L,y,N,N') = L^P U (N,N'). (469) 
B. Action of the Poincare group on vacuum 

Vacuum states are neither exactly translation invariant (in consequence of the "zero 
energy" 4-momentum, which we do not neglect) 

U(l,y,N,N')\0,N,N') 
= e^ N ^\0,N,N') 

= ( J dpe iyp O m (p)\p) ® |O,O)) 0JV ® ( J dke iyk/2 O {k)\k) ® |0) X ^ 
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nor Lorentz invariant 



U(L, 0, N, N')\0, N, N') = (u m (L, 0, lf N ® U (L, 0, lf N ') (\O m , 1)* N ® \O , lf N ') 

= [C/ m (L,0,l)|O m ,l)J ® (U (L,0,l)\O ,l)) 
C/ m (L,0,l)|O m ,l) = y dpIpX^-Vl^Cl®!) y dp'OmCpOIP'J^lO.O) 



y dp'o m (p / )ip)(^"W>«'io,o) 

dpOrniL-^p)® |0,0), 

C/ (L,0,l)|O ,l) = y dfcjfcXi- 1 *;! <g) 1 y dk' O (k')\k') <g> |0) 
= J dk J dk' O (k')\p) (L^klk') <g> |0) 
dA;Oo(^ _1 fe)|fe) <g> |0). 



We can say that vacuum wave functions behave under Poincare transformations as classical 
scalar fields on mass-m hyperboloid and light cone, respectively, 

O m (p) e^O m (p), (470) 

O m (p) O m (X- 1 p), (471) 

Oo(fe) O (i- 1 fe), (472) 

Oo(fe) -)• e^ fe/2 O (fc). (473) 

The "zero-energy" parts of 4-momentum can be removed by a well defined unitary trans- 
formation. Such a new vacuum will become 4-translation invariant, but one cannot do the 
same with Lorentz transformations. Lorentz invariance would require constant O m (p) and 
O (k), a condition inconsistent with square-integrability 

y dp\O m (p)\ 2 = y dk\O (k)\ 2 = 1 (474) 

which we assume. 

The entire subspace of vacuum states is nevertheless Poincare invariant. It follows that 
the projector on the vacuum subspace 

/ r \®n , r x ®n' 

U (N,N') = U dp|p)(p|(g)|0,0)(0,0|J ® U dk\k){k\ ® |0)(0|J (475) 
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is Poincare invariant 



U(L,y,N,Nyn (N,N r )U{L,y,N,N') = U (N,N'), (476) 
and commutes with U (L, y, N, N'), 

U (N,N')U(L,y,N,N') = U(L, y, N, N')U {N, N'). (477) 

C. Vacuum in 4-position space and nontrivial boundary conditions 

The fact that vacuum wave functions behave under Lorentz transformations as classical 
scalar fields on mass-m hyperboloid and light cone, respectively, suggests their interpretation 
as Fourier-space amplitudes of classical scalar fields 

O m (x) = j dp(0Me- vpx + 0jp)e^, (478) 
O {x) = J dk (O {k)e~ ikx/2 + 7Mkje lkx/2 }. (479) 

The fields Oq(x) and O m (x) provide a means of imposing nontrivial boundary conditions 
associated with concrete geometry. 

This is an interesting point that will be discussed later in the context of the Casimir effect 
and fields in cavities. At this place it is just good to know that boundary conditions should 
be in my approach imposed on both vacuum states and field operators. 

Remark: The essential mathematical difference between \0,N,N') and the vacuum |0) of 
the bosonic Fock space is related to Poincare invariance. Uniqueness and Poincare invari- 
ance of the vacuum state is one of those Wightman axioms of standard Fock-space based 



quantum field theory [45|, |46J that are not satisfied by my formalism. (Note, however, that 
vacuum understood as the subspace of vacuum states is both unique and Poincare invariant 
even in my formalism; in Wightmanian formalism this space is 1-dimensional so both types 
of invariance and uniqueness are then equivalent.) In the context of indefinite-frequency 
oscillators the state \0,N, N') represents a gas of N massive and N' massless bosonic oscil- 
lator wave packets at zero temperature. As such it should be neither unique nor Poincare 
invariant, similarly to states of actual Bose-Einstein condensates occurring in atomic physics 
(clouds of trapped cold atoms are not even translation invariant). In relativistic theories all 
states should be at least Poincare covariant.k 
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Remark: In solid-state physics one distinguishes between "particles", such as atoms forming 
a given medium, and "quasi-particles" (phonons, plasmons, etc.) - - quantized oscillations 
of the medium. The physical difference between the Bose-Einstein condensate \0,N,N') 
and the Fock vacuum |0) can be better understood if one thinks of Lewenstein-You-Cooper- 
Burnett Fock-space formulation of atomic Bose-Einstein condensates |47J]. The vacuum there 
corresponds to the Fockian |0), i.e. the state of no atoms. This is a purely formal object, 
a "ground state" meaning "no atoms in the trap". A "Bose-Einstein condensate at zero 
temperature" is a state analogous to \0, N, N'). It is not necessarily the lowest energy state 
of the total Hamiltonian since electrons in atoms may exist in higher excited states, but 
rather a ground state of a part of the Hamiltonian (typically representing the center-of- 
mass kinetic energy). In my formalism a kind of |0) could correspond to N = 0, N' = 
(in the first version of the formalism published in |7J I indeed considered this type of dual 
particle/quasi-particle Fock-type structure), n-quasi-particle states, that is those containing 
n excitations of the oscillators, are regarded in these notes as representing n particles. The 
notion of particle statistics refers only to the quasi-particle level. Numbers of quasi-particles 
may change in time as opposed to iV and N' that remain fixed. The "medium" composed of 
our N and N' particles is static and does not explicitly take part in dynamical processes. A 

Remark: The so called loop diagrams I will describe in detail later correspond to "vacuum 
to vacuum" processes, such as annihilation of particles that were spontaneously created 
from vacuum. In my formalism the appropriate formulas will be obtained by sandwiching 
evolution operators between projectors on vacuum and not between vacuum state vectors, 
as would be the case in more standard approaches. A 

XVI. INTERACTION-PICTURE DYNAMICS 

Let (f) be any operator and H = Hq + Hi a Hamiltonian. Its dynamics in Heisenberg 
picture is given by 

^(t) = e iHt <f)e- im = U 1 (ty e iHot <J)e- iHot U 1 (t), 0(0) = 0, (480) 
U x {t) = e iHot e- im , (481) 
iJJ^t) = e iHot {-H + H)e~ im = e^H^-^U^t). (482) 
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Denoting 

(j)(t) = e lHot <j)e- lHot , (483) 

Hi(t) = e iHot H ie - iHo \ (484) 

we can write 

4>(t) = U^^U^t), (485) 

itliit) = H^U^t). (486) 

Let us note that <f>(t) and <f>(t) are equal at t — — this is an example of initial condition 
for Heisenberg dynamics. It is often convenient to take the initial condition 0(to) = (f>(to) at 
an arbitrary t , 

4>(t) = e iHit - to U(to)e- iHit - to) (487) 

= e iH( ~ t ~ t ^ e iH o t o ( j )e -' iH o t o e - iH ( t - t o) (488) 

_ e iH{t-t ) e -iH (t-to)^^ e iH (t-t ) e -iH{t-t ) (489) 

= Vi(Mo)V(W(Mo). (490) 



The evolution operator 



V 1 {t,t Q ) = e iHo(t-to) e -iH{t-t )^ ( 491 ) 



satisfies Schrodinger's equation 
at 

= e iH (t-to) jj -iH(t-t ) 

= e iHo(t-to) Hie -iH (t-to) Vi ( t ^ tQ j 

= ifi(t-* )^i(Mo), (492) 

and 

Vi(to,*o) = I (493) 
Vi(Mo) = ^(*-* ,0), (494) 

but in general does not fulfill the composition property, i.e. 

^i(ti,t 2 )V 1 (t 2 ,t 3 ) ^ V 1 (t 1 ,t 3 ). (495) 
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Another option is to take 

Ux(t,t ) = jHot e -iH(t-t ) e -iH t ^gg) 

which, as opposed to Vi(t,to), satisfies the composition property 

C/i(ti,t 2 )t/i(t 2 ,t 3 ) = «/i(ti,t 3 ). (497) 
The dynamics is now different: 

i^-Ui{t,to) = - e iHot Hoe -iH(t-t ) e -iH t + e iH t He -iH(t-t ) e -iH t 
CI I 

= HiWUfato). (498) 
The initial condition is unchanged, 

£M*o,*o) = I, (499) 
but Ui(t, to) 7^ U\(t — to, 0). Assuming 0(io) = we find 

0(t) = e i^(*-*o)^ e -iff(t-to) 

= e-^Utit, toYe^fe-^U^t, t )e iHoto 

= e- iHo ' f/ 1 (t,t )V(t)f/i(t,to)e iHo * (500) 

Let us note that 4>(t ) = <fi = 0(0) whereas the case of Vi(t, t ) corresponded to 4>(t ) = (fi(t ). 
For to = the three evolution operators, (I48ip . (I49ip . (I496p . coincide: 

U 1 (t) = U 1 (t,0) = V 1 (t,0). (501) 

Equation (I498P is known as the interaction or Dirac picture. (149 8 p and (1499 p are equivalent 
to 

E/i(t,t ) = I+H) / dTi#i(ri)EMri,*o) 

= i + H) / dn^^rOfi+H) f 1 d^H^U^to)). 

Jtn V -/io 7 
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By iteration 

l7i(Mo) = I+H) / rfr 1 iJ 1 (r 1 ) + (-2) 2 / dn / 1 rfr 2 iJ 1 (r 1 )iJ 1 (r 2 )+... (502) 

•/to -'to -/to v v ' 

loop terms originate here 

= Texp ( - % J dn #i(ti)) . (503) 

Texp(. . . ) is the so-called time-ordered exponential. 

Another popular (and equivalent) form of Texp(. . . ) involves the so called time-ordered 
products, 

tA(Mo) = jt(-*T [ dri j 1 dr 2 ... r 1 dr n H 1 (T 1 )H 1 (r 2 )...H 1 (T n ) 

00 nt ft 

= Z)H) n / dn... dr n 9(T 1 -T 2 )...9(r n ^-r n )H 1 (T 1 )...H 1 (r n ).(50A) 

Now we have n integrals from t to t. Taking into account all the possible permutations 
(Tj 15 . . . , Tj n ) of the n-tuple (ri, . . . , r n ) and assuming that orders of integration can be inter- 
changed (which is not always obvious, especially if to or t are taken at infinity) we obtain 

l7i(Mo) = f^^T- I dn-.. f dr n T(H l {r 1 )...H 1 {r n )), (505) 

n=0 n ' ^ -/to V 7 

where the time-ordered product is defined by 

T(H l {r 1 )...H l {r n )) = £ 0(r n - r l2 ) . . . 9^ - njH^) . . . H^nJ. (50Q) 

The above forms of dynamics can be further generalized as follows. Assume O C 1R 4 is 
a subset of Minkowski space contained between two spacelike hyper-surfaces E n and £ T2 . 
Let us also assume that there exists a 1-parameter family of hypersurfaces £ r , O = [j T E T , 
such that for each x G there exists one and only one E T containing x. The family S T is 
termed the foliation of O. 

Interaction-picture Hamiltonian associated with a given foliation is defined by 



V(t) = / d£ T (x)V(x). (507) 
V(x) is a Poincare covariant interaction-Hamiltonian density, 

U(L, y yV(x)U(L,y) = V{L-\x-y)). (508) 
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U(L,y) is a representation of the Poincare group — at this moment we do not have to 
specify exactly which representation we have in mind. dH T is a measure on S T satisfying 

d A x = drdi: T . (509) 

Condition f !509j) implicitly enforces relativistic covariance of the formalism. Let x' = Lx 
where L is a Lorentz transformation. Then 

d 4 x = dxa d 3 x (510) 

= dx' d 3 x' (511) 

= dr - ^ = dr - d3x ' (512) 
^/1 + x 2 /t 2 ^1 + x' 2 /t 2 

define two hyperplane foliations E r = {x; xq = r} and E T = {x; x' = r} of the Minkowski 
space, and the hyperbolic foliation S r = {x; (x — y) 2 = r 2 , x > y } of the future causal 



cone of some event y (the so-called Milne universe 48j). The respective measures are 



dH T (x) = d 3 x, (513) 

dE T (x') = d 3 x', (514) 

aEJx) = — = = — = = dzlJx ). (515) 

^1 + x 2 /t 2 y/l + x' 2 /r 2 



In Dirac's terminology 49] the first two foliations define instant-form dynamics. The hyper- 
bolic foliation defines the point-form dynamics in the Milne universe. Point-form quantum 
optics of classical sources formulated in "my" approach to quantum field theory was de- 



scribed in detail in 



Milne universe is the exceptional case of dynamics where free-field 
initial condition at r = is physical: At r = the fields are at the boundary of the universe. 
Interaction-picture evolution operator with respect to any foliation is a solution of 

i^-u x {T,n) = v^u^n), (516) 

dr 

U^run) = I. (517) 

An important property of interaction-picture dynamics is the fact that Hamiltonians H\(t) 
are constructed from operators that depend on time through the free-Hamiltonian dynamics 
generated by H . 
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XVII. CONCRETE CHOICE OF V 



A scalar-scalar interaction that seems formally closest to quantum electrodynamics is 

V{r) = I d3X 2/ -j»(x,N,N')<j>(x,N,N>) 
Jr3 y/i + x 2 /r 2 r 

d x 

x»jp(x,N,N')<l>(x,N,N') (518) 



x z 



with the current given by 

j„(x, N, N') = iqip(x } N, N'^d^ix, N, N') - iqd^ix, N, N')^(x, N, N'). (519) 



The point x = (V r 2 + x 2 , x) belongs to S T = {x; x 2 = r 2 , xo > 0}, so 

x/t = ( v 7 ! + x 2 /t 2 , x/t) -> (1, 0) withr^oo. (520) 

As our first nontrivial exercise, leaving aside for the moment the issues of relativistic covari- 
ance but concentrating entirely on the problem of divergences in loop diagrams, we begin 
with a simplified version of (1518j) . 

V(t) = [ d z x J{i {T,x,N,N')<P{T,x,N,N'). (521) 



It corresponds to the asymptotic form of ( 15181) for r — > oo. 

One should not confuse the models with the so-called scalar electrodynamics, a gauge 
invariant theory describing interaction of charged spin-0 particles with electromagnetic fields. 
The corresponding interaction Hamiltonian is much more complicated but contains a term 
analogous to ( I518p (see example 8.6 on p. 251 in |50|). 

XVIII. SIMPLIFIED HAMILTONIAN (l52T|) IN MOMENTUM SPACE 

We begin with expressing (152 ip in momentum space. It is this place where we have to 
carefully define products and integrals of operators [see the discussion accompanying (I236p 
and f !353p ]. The integrals involve dp = d 3 p/[(2Tr) 3 2^m 2 + p 2 ], dr = d 3 r/[(2n) 3 2y/m 2 + r 2 ], 
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dk = d 3 A;/[(27r) 3 2|fc|], p = ^m 2 +p 2 , r = Vm 2 + r 2 , k = \k\. 

V(r) = q J f X J A P {a {p )e-»* + « P )V) iBo J + 

-* j *{<***-<» + •wv)/#(wr- + aW ^>w 

= 1 J ^ X J dpdrdkr 

— g y d 3 a; ^ dpdrdkpo 

x (a(p)e-* p3: - &(p) t e ip!e ) ^(r)e" ira + a(r) t e ira ) (c(A;) e - fc + c^e***). 

The above expressions should be understood in the following sense 

V(t) = q J d 3 x J dpdrdkr Q a(p)b(r)c{k)e- l{p+r+k)x + . . . 

= lim q / d 3 x / dpdrdkr® 

ni,n2,ri3->-oo J J 

xa(p, m, iV, JV')&(r, n 2 , iV, JV')c(A;, n 3 , JV, jV')e- i(p+r+fc)ie + . . . 
Here a(p, ni, A/", A/"') is the (N, AT^-oscillator extension of 

a(p,ni,l) = \p, n 1 )(p, ni| ® (a <g) 1) (522) 



or 



a(p, ni,^,!) = |p, ni)(p, n'J <g> (a (g) 1) (523) 
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(compare the discussion of spectral theorem for M-shaped Dirac deltas). Similar definitions 
apply to the other operators. For finite rii, n 2) ... we can apply (I212p so that 

V(t) = | J ^ 3 rrfA;a(p)6(r)c(fc)e~ l( v /m2 +P 2 +^ m2+ ^+l fc l) T ( 5( 3 )(-r- - p - k) 
-\ \ rfV^^a(p)^(^)c(fc)e~ l( v /m2+p2+Vm ' +r2+ l fc l) T 5( 3 )(-p -r-k) 
+ | / ^ 3 rdfca(p)6(r)c(fc) t e- i( v /m2 +P 2 +^+^-|fc|)^(3)(_ r _ p+ fe ) 

- - / rfVrrfA;a(p)6(r)c(fc) t e- i (v /m2 +P 2 +^+^"|fc|)r (5 (3) ( _ p _ p + fc) 

-|/« o(p)oW , (fe)e -«™,^ ( _, +p+fe) 

- * /Wa( P )a(r)tc( fe )e-'«-v^ +W ^< 3 . ( _ p + r _ fc) 

+ | / ^VrfA;6(p) t 6(r)c(fc)e- i( -v /m2 +P 2 +v /m2 +^+l fc l)^( 3 )(_ r + p _ fe) 
+ | y rfVrrfA;6(p) t 6(r)c(fc)e- i( -V^ 2 +P 2 +^+^+|fe|)r 5 (3)(_ p + r + fc ) 

+ | / ^ 3 r^6(p) t 6(r)c(fc) t e- i( -v /m2 +P 2 +^+^-l fc D^(3)(- r + p + fc) 

+ §/ ^V^Kp) t ^(^)c(fc) t e^v /m2 +P 2 +^ /m2 +^-l fc l)^( 3 )(-p + r- - fc) 

-|/« wtaWtc(fe)e -^^ +l ,^ ( _ r _ p+fc) 

+ | / rfVr^6(p) t a(r) t c(fc)e- j( -v / - 2 +P 2 -^+^+|fc|)r 5 (3)(_ p _ r + fc ) 

- ^ ^rf 3 rrfA;6(p) t a(r) t c(fc) t e- i( -v /m2 +P 2 -V m2 +- 2 -l fe l)^( 3 )(-r -p- k) 
+ | / d 3 pdr^6(p) t a(r) t c(fe) t e- i( -V" l2 +P 2 -^+^-|fc|)T 5 (3)(_ p _ r _ fc) 
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V(r) 



+ 
+ 



+ 1 



+ 
+ 



+ 1 



+ 



Q 



dpdk (a(p)b(-p - fe) - a(-p - fc)6(p))c(fc)e- i( v /m2 +P 2 +v /m2 +(p+ fc ) 2 +l fc l)- 
dpdk (b(-p - fc) t a(p) t - b( P y a (-p - fc)t) c (fc)t e l (v / ™ 2 +P 2 +v /m2 +(p+ fc ) 2 +l fe l)- 
dpdk (a(p)b(-p + fe) - a(-p + fc)6(p))c(fc) t e- i( v /m2 +P 2 +v /m2 +(p- fc ) 2 -l fc D- 
dpdA; (&(-p + fc) t a(p) t - &(p) f a(-p + fc) t )c(fc)e i( v /m2 +P 2 +v /m2 +(p- fc ) 2 -l fc D- 
dpdA; (b(p + fc) f 6(p) - a(p)a(p + fc)t ) c (fc) e -*(V m2 +P 2 V m2 +(p+ fc ) 2 +l fc l)^ 
dpdA; ffe(p) t 6(p+ fe) - a(p + fc)a(p) t )c(fc) t e t( v /m2+p2 -v /m2+(p+fc)2+|fc|)r 
dpdk (b(p) ] b(p - fe) - a(p - fc)g(p) t )c(fc)e- t (-v /m2+p2+ v /m2+ ( p - fc ) 2+ l fc l) r 
dpdA; (&(p - fc) f 6(p) - a(p)a(p - fc) t )c(fc)¥ ( -v' m2+p2+ v / " 



/m 2 + (p-fc) 2 + |fc|)r 



I have underlined the four terms that are not normally ordered. Defining the normally 
ordered part of V(r) by 



:V(r) 



1/ 



dpd/c 



+ 2 / 



+ 



+ 



1/ 



dpd/c 
dpd/c 



+ 



+ 



+ 2 / ^P 6 ^ 



1/ 
/■ 



dpd/c 
dpd/c 



- ^ / dpd/c 



a(p)b(-p - fe) - a(-p - fc)6(p))c(fe)e- l (v / ™ 2 +P 2 +v /m2 +(p+ fc ) 2 +l fe l)- 
6(-p - fc) t a(p) t - 6(p) f a(-p - fe) t )c(fc) t e i( v /m2 +P 2 +v /m2 +(p+ fc ) 2 +l fc l) 7 



a(p)fe(-p + fe) - a(-p + fe)6(p))c(fc) t e- l( v /m2 +P 2 +v /m2 +(p- fc ) 2 -l fe l)- 
6(-p + fc) t a(p) t - 6(p) t a(-p + fe) t )c(fc)e l (v / ™ 2 +P 2 +v /m2 +(f- fc ) 2 -l fc l)- 
6(p + fc) t 6(p) - a(p + fc)t a (p) ) c (fc) e -^(v / ^ 2 +p 2 -v / " 2 +(p+fc) 2 +|felK 
6(p) f 6(p + fe) - a( P y a (p + fc) N )c(fc) t e^v / ' m2 +P 2 -v /m2 +(p+ fc ) 2 +l fc l)- 



b(p)%p - fe) - a(p)tq(p - fc) )c(fe)e- t (-v /m2+p2+ v /m2+ ( p - fc ) 2+ l fc l) r 
b(p - fc) t 6(p) - o(p--fc)ta(p)) c(fc)t e i (-v /m2 +P 2 +v / 



m 2 +(p-fc) 2 + |fc|)r 
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where the underlined terms are reordered so that annihilation operators are to the right of 
creation operators, one finds 

V(t) = : V(t) : 

-| J dpdkI m (p)2(2n) 3 v ^ 2 Tp^S* {3 \k)c(k)e^^ ; ^ 

-| J dpdkI m (p)2(2Ttf^m 2 + p 2 ( j*(3)( fe ) c ( fe )t e i(V" l2 +P 2 -V" l2 +Cp+ fc ) 2 +l fc l)^ 
-| y dpdkl m (p)2(2n) 3 v / m 2 + p ^*( 3 )(k)c(k)e-^-V m2 +P 2 W m2 +(p- k ) 2 +\ k \> 
-| y ^A;J m (p)2(27r)V^ 2 + P 2 ^ (3 Hfc)<fc) t e i( "^ m2+p2+v/m2+(p_fc)2+|fc|)T 



V(r) 



(524) 

Let us note that in the standard approach, where I m (p) = 1, the vacuum term 

V(r)- : V(r) : (525) 



/ — -— ^==/ m (p)(2 7 r) 3 2 v / ^T^= / rfV (526) 



is badly divergent due to the first integral 

d 3 p 

(2n) 3 2^m 2 + p 2 

This is precisely an example of "ultraviolet catastrophe" since the integral is divergent due 
to the behavior of its integrand (equal to 1) for large momenta. 
The other term 

/ OT i<3,(fe) ( c(fe)+ ' :(fe)t ) ^Slw (527) 

would be in standard approaches an example of "infrared catastrophe". Eq. f)527p shows 
that c(k) + c(ky must tend to zero at least as fast as |fe| but in the standard formalism 
it is completely unclear how to achieve it in a mathematically precise way. The "infrared 
catastrophe" is typical of massless quantum fields — |fc| in denominators would be replaced 
for m > by \/m 2 + k 2 . 

The most typical way of dealing with the divergences occurring in V(r)— : V(r) : is just 
to ignore them! We ignore them not because they are small, but "because we do not want 
them". One replaces V(r) by : V(r) : and starts with 

i^pU^rx) = : V(t) : Ux(t,Tx). (528) 
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It is not hard to guess that this type of "solution" will sooner or later lead to another 
"catastrophe" that will have to be remedied in one way or another. 

The elegance typical of quantum mechanics has been completely lost. What one does is 
not "sensible mathematics" . This is not even a physical theory but a "collection of working 
rules" — at least in Dirac's opinion. 

Let us now see what happens in (JV, N') representations of HOLA. First of all, by (140 2 j) 



I m (p,N,N') = I m (p,N)®I {N') 
1 

~ N 

The expression 



Tf((\p)(p\® (1® 1)) ®/ m (l)® J TO (1) + ...) ®I Q (N'), 



I 7Z- -^===UP,N)(2nf2V^Tp~ 2 (529) 
J {Zny ; 2y m 2 + p 2 

is l/N times a sum of N operators of the form 

/ - - J m (p,l)(27r) 3 2v/m 2 + p 2 = (2vr) 3 2 / dp \p)(p\ ^m 2 + p 2 <g, (1 ® 1) 

J \2tx) a 2^/ m 1 + p 2 J 

(530) 

each of them acting in the Hilbert space of a single indefinite-frequency oscillator. (15301) is 
the spectral representation of a well behaved operator 



(27r) 3 2 v / m 2 +p 2 g> (1 g> 1), (531) 

where 

p = J dp\p)(p\p. (532) 

Recall that (p|p) = 0, a fact that will become important later. 

The example shows "regularization by quantization" in action: There is no "ultraviolet 
catastrophe" because the ordinary integral 

jdp{...) (533) 

has been replaced by the spectral integral 

dp\p)(p\(...) (534) 
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Operators (153 ip and (I532j) are the analogues of, respectively, 

f2 = u)\u) (cu\ (g> /, (compare (1971) ) 

and Cj = 2~2ui UJ \ UJ ) ( u \ we nave started with in our discussion of indefinite-frequency oscilla- 
tors. 

Now let us discuss the issue of infrared divergence of 

d 3 k 



(27r) 3 2|fc 

d 3 k 



<f {3) (k) (c{k, N, N') + c{k, N, iV / ) t ' 
UN) ® / f2 ^ |fc| ^ (3) (fc)(c(fc,iVO+c(fc,iVO t ). (535) 



(27r) 3 2|fc 

This is essentially a sum of N' single-oscillator operators 

/(27^^ (3)(fc) ( c(fc ' 1) + c(fc ' 1)t ) = /(^||^^ 3) ( fc )l fc )^l®( G + at )( 536 ) 

multiplied by Now let us act with it on a single-oscillator (N' = 1) vacuum state 

d 3 k 



(0 k .S*Q\k)\k)(k\®(a + <f)\O ,l) = J dk5* {3 \k)\k)(k\®a^ J dk' O {k')\k', 0) 

5^ 3 \k)O (k)\k,l) 



(2vr) 3 2|fc 



= lim 9 1 O (fc)|fc,l) = 0. 
fc->o {2ir) 6 2\k\ 

The latter follows from the Poincare invariant boundary condition (144 ip 

fc^O |fe| n 

Now check the action of (I536p on a more general "one-photon" one-oscillator state 

= J dk'^(k')c(k',iy\Oo,l), (537) 



d 3 k 



+ 



^|^5*( 3 )(fc)c(fc, l)t | dk^(k')c(k', l)t|O , 1) 

5* (3) (fc) / dk'^(k')5 (k, k')I (k, 1) / dk"O {k")\k") <g> |0> 



/ 

+ | j2^\k\ S<3){k)c{kA)t J MMWrfJ dk"O (k")\k")®\0) 
/ (2 ^32|fc| ^ (3)(fc) / dk '^ k 'Mk,k') J dk"O (k")\k)(k\k")®\0) 
+ | (2 ^ |fc| ^ (3) (fc) / ^V(fc') | rfA;"Oo(fc // )|fc)(fc|fc'}(fc'|fc")®|2) 
(270 3 2|fe| ^ (3)(fe) / dk '^ k ') 5 ^ k ') j dk"O (k")5 (k,k")\k)®\0) 
/ (270 3 2|fc| <r(3)(fc) / / rfA;//0 o( fc/, )^( fc ' fc/ )^( fc/ ' fc// )l fc )®l 2 ) 



+ 

= / (2^q^ (3) (fcMfc)Oo(fc)|fc)®|0) + | ^|^^ 3 )(fcMfc)O (fc)|fc)®|2) 

= L m o(27)W^ )0o(fc)|fc ^( |0) + |2) )- 

The term again vanishes unless blows up at the origin faster than any power of l/\k\. 
The argument can be generalized by induction to states generated from vacuum by any 
number of creation oscillators. 

All of this is possible since in the formula 

O (k) = (O ,N'\I (k,N')\O ,N') (538) 

the operator Io(k, N') is not proportional to the identity. In representations occurring in the 
standard formalisms one would find similar expressions as we have obtained (since algebraic 
manipulations with operators would be identical), but with 

Oo(fc) = (0\I(k)\0) = (0\I\0) = 1 (539) 

and then infrared divergences would be inevitable. 
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FIG. 6: The simplest loop vacuum-to- vacuum Feynman diagram represents a process where cre- 
ation operators a{p)\ b(—p — fc)t, c(fe)t act on vacuum and the resulting state is then acted upon 
by the annihilation operators a(p'), b(—p' — k'), c(fc'). So one starts from vacuum and ends up 
with vacuum. 

XIX. VACUUM-TO- VACUUM LOOP DIAGRAM IN 2ND ORDER PERTURBA- 
TION THEORY 

Vacuum-to-vacuum loop diagram is the 2nd order perturbative correction to 
(O, N, N'\Ui(t2, ti)\0, N, N'}, where \0,N,N'} is a vacuum space. What is very impor- 
tant we want to have a result that does not depend on concrete explicit forms of vacuum 
wave functions O m (p) and Oo(fe) but only on their general properties such as boundary 
conditions. 

So let n = n (iV, N') be the projector on the vacuum subspace. It is characterized by 

a(p)U = b(p)U = c(fc)n = n a(p)t = U b( P y = n c(fc)t = 0, (540) 
J dk5* (3) (k)c(k)^U = J dk5* (3) (k)U c(k) = J dk5* {3 \k)I (k)U = 0. (541) 

Our goal is to compute in 2nd order perturbation theory the operator iloL r i(r 2 , Ti)n . The 
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appropriate term is 

U V(r)V(r')U 

= n : V(t) :: V(r') : n 

= — J dpdkdp' dk' e^^^^+V™^ 

x n (a(p)&(— p — k) — a(—p — k)b(p)jc(k) 

x (b(-p' - fc / ) t a(p / ) t - b( P ya(-p' - kyy(kyu 

2 p 

= — dpdkdp e - i (V m2 +P 2 +V m2 +(p+ fc ) 2 +l fc l) T e i( V m2+p ' 2+ V m2+ (p'+ fc ) 2+ l fe l) r ' 

x U Q (a(p)a{p')%-p - k)b(-p' - fc) f + a{-p - k)a{-p' - k) ] b(p)b(p') ] 

- a{p)a{-p' - k) ] b(-p - k)b( P y - a(-p - k)a(p') ] b(p)b(-p' - fc) t )/ (fc)n 

= — dpdkdp e - i (V m2 +P 2 +v /m2 +(p+ fc ) 2 +l fc l) T e i (v /m2+p ' 2+ v /m2 +(p'+ fe ) 2 +l fe l) T ' 

+ I m (-p - k)5 m (-p - k, -p - k)I m (p)5 m (p,p') 

- I m (p)S m (p, -p' - k)I m (-p - k)S m (-p - k,p') 

- I m (-p - k)5 m (-p - k,p')I m (p)5 m (p, -p - k) ) J (fe)n (542) 



The latter expression is ill denned in standard approaches due to occurrences of squared 
Dirac deltas. Additional divergences would be generated if one assumed the standard condi- 
tion Io{k) = I m {p) = 1. This is why the loop diagram is normally an example of a "difficult 
infinity" 
4.1 of 



see for example an analogous calculation made in full electrodynamics in chapter 
5 ljl). In my formalism the calculation is rather trivial. 
We have to keep in mind that dk = p (k)d 3 k, dp = p m (p)d 3 p so that internal consistency 
conditions analogous to (I284p can be derived. The only solution is thus to work with M- 
shaped Dirac deltas satisfying 6 m (p, p) = 0. The result is 

x [5 m (-p- k, -p-k) + S m (-p - k,-p- k) - 5 m (p,p) - 5 m (p,p)^j 
x U I m (p)I m (-p - fe)J (fe)n Q 
= 0. (543) 

Notice that we have not needed the explicit forms of O m {p) and O (k). In standard quantum 
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field theoretic parlance the result thus "does not depend on a cutoff' . 



XX. MASSLESS QUANTUM SCALAR FIELD PRODUCED BY A CLASSICAL 
POINTLIKE SOURCE 

In order to compute radiative corrections to atomic energy levels one first has to decide 
what kind of potential is associated with atomic nucleus. In simplest approaches the nucleus 
is modeled by the classical static pointlike charge density 

p{x) = q5 {3) (x). (544) 

Classically the field produced by the source would be a solution of an inhomogeneous 
(Maxwell, d'Alembert...) field equation. In the quantum context one should solve the 
Heisenberg equation with appropriate interaction term. 

Let us employ Heisenberg and Dirac picture forms of evolution equations to the dynamics 
of 4>(t,x) for V(t) = f R3 d 3 x p(x)(()(t,x), 

<j>(t,x) = J dA;(c(fc)e- < l fc l* +ifc - a, + c(Jb) t e i l fc l t - ifc - SB ), (545) 

V(t) = q [ d 3 x5 (3 \x)(f)(t,x) 

= #(*,0) (546) 

= q J dkc(k)e' ilklt + q J dk c(fc) f e j|fc| * . (547) 

s v ' v v ' 

V-it) V+(t)=V_(t)t 

The example is interesting since in standard field quantization Ui(t,t ) does not exist due 
to ultraviolet divergences. We will need 

[V-(t),V + (t')} = q 2 [J dkcitye-'W, J dk'cik'^eW' 

= q 2 j rfA;/ (fc)e i|fc|(t '~ t) . (548) 

Of particular interest is the equal-time commutator 

[V_(t),V + (t)\ = q 2 J dkl (k). 

In reducible iV-representation of HOLA 

[V-(t),V + (t)] = q 2 J dkI (k,N) (549) 
= g 2 / (AO (550) 
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is a well defined operator (J549]). This should be contrasted with the divergent integral that 
would have occurred here for Jo(fe) = 1. Since for all times 

[V4t),V„(t')] = [V + (t),V + (t')} = (551) 

and 

[\V-{t), V + (t%V ± (f)] = 0, (552) 

we can use the continuous Baker- Campbell-Hausdorff (BCH) formula (cf. Appendix H, 



Eq. (13) in |52|), 



l/i(Mo) = Texp( -ijf dr(F + (r)+V_(r))) (553) 

= exp ^ — z J dr V+(t) j exp ^ — z y dr K_(r)^ 

xexp( / dr / dr'fK^r'X'Mr)]). (554) 

Another useful form of continuous BCH identity is based on the ordinary BCH identity 
e x e Y = e lX,Y ^ 2 e x+Y ^ which can be applied to the right-hand-side of (15541) . 

EM*,*o) = exp(-z f dTV(r))exp(-l- f dr f dr' ' [V+(r') ,F_(r)]) 

V J to / \ I J to J tQ / 

xexpt [ dr [ dr' 6{r - r')[V + {r'), VL(r)]) (555) 
= exp(-i f drV(r))exp(- [ dr [ dr' sgn (r - r')[V + (r'), VL(r)]) . 



(556) 



sgn(x) and 6(x) are, respectively, the sign-of-x and step functions. For r = r' the commuta- 
tors in exponents are finite in "my" representation of HOLA, so exact values at zero of sgn 
and 9 are for the moment irrelevant. Both exponents at the right side of f )556D are unitary. 
The second of them commutes with all the elements of HOLA, so this is an operator that 
in practice behaves as numerical phase factor. 

Continuous BCH formula can be verified in a straightforward manner. First of all, at 
t = to the left and right sides of (I554p equal I. It remains to check if they satisfy the same 
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differential equation, 
.d£/i(Mo) 



dt 



V+(t) exp ( - i jf dr V+(r) ) exp ( - i jf" dr V_(r 
xexp( /" dr [ T dr'[V + (r'),V^r)]) 

V -/to 7 

+ exp(-i / dr V+(r)V_(i)exp ( - % I drV-{ 
xexp( / dr f dT'[V + {T'),V^T)]) 

+i [ dT , \y+{T , ),V-(t)]exp(-i [ drV+(r)) exp (-i [ drV_(r)) 

Jt ^ Jto Jt ' 

xexp(/ dr f dr'[V + (T'),V^T)} 

V Jto Jto 



V+(t) + e 



to J to 

iSt drV + {r) 



V-(t)e if !o dTV+(T) +i / dr'[V + (r'),V4t)] }U 1 (t,t ) 

to 



Employing ( 1531) and ( 1552R we rewrite the second term 



iJl<iTV + (T) 



e Jt o 



if! drV+(r) 



VC(t) + -i \ drV+(r),V4t) 

Jto 



which ends the proof: 



We can finally compute 



iUi(t,t Q ) = (V+ty + V-ityU^to). 



(557) 



4>(t,x) = exp [i I drV(r) )(j)(t, x) exp 



to 



t 

i, / drV(r] 

to 

t 



exp (iq J dr 0(r, O)^0(i, a;) exp ^ — iq J dr0(r, 0) 



'to 

4>(t,x)+iq dr[(j)(T,0),(j)(t,x)] 
Jt 

<f>(t,x) 
ft 



+iq J dr J dkdk' 

x [c(fc)e- i|fc|r + c(fc)V |fe|r , c(fc')e- i|fc ' l ' +ifc '- a; + c(fc') t e i| *' l *- < * ,,fls ] 
0(t,ai)+ig / dfc/ (fc) (e^ 1 *"** / dre~ mT - e~ mt+ ' lk - x [ dre 

J Jt Jt 

<f>(t,x) 



i\k\r 



-Q 



d 3 k 



(2vr) 3 2|A;| 



(j){t,x) - q 



Mk)(e 

d 3 k 



(2vr) 3 2|fc| 



i|fe|t— ik x / — i\k\t — i|fe|to 



) + 



i\k\t+ik-x ^ p i\k\t p i|fe|to 



-Io(k) ( e~ ik ' x -\- e ik x _ e 'ikx e -i\k\t _ e -ikx e i\k\to 
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(f)(t, x) splits into free-field and source parts. The source term can be further split into the 
"generalized Coulomb field" 

/d 3 k 
(27r) 3 \k\ 2 l °^ COSk ' X ^ 558 ^ 

and the "compensating field" 

/d 3 k 
^ 3 | fc | 2 A)(fc) cos(A;a; - |fc|* ). (559) 

Compensating and Coulomb fields cancel each another at t — t . Let us discuss the two 
terms separately. 

A. Compensating field 

In order to better understand the role of the compensating field let us consider an arbitrary 
coherent state 

\a,N) = V (a,N)\O ,N), (560) 
V (a,N) = exp J dk(a(k)c(k, iV) f -a(k)c(k, N)j. (561) 

Since 

V (a,N^I (k,N)V (a,N) = I (k,N) (562) 

one finds 

/d 3 k 
(27r)3|fc|2 (Qo, N\I (k, N)\O , N) cos(kx - \k\t ) 

= q j (2^f\k\^ ° Q<yk ^ 2 C ° S<ykX ~ \ k \ t °)' 

The same result will be obtained if we replace the displacement operator by any unitary oper- 
ator constructed from elements of HOLA. Normalization and boundary conditions imposed 
on vacuum imply 

/ v$m mk)? = 11 (563) 

lim^ = 0. (564) 

fc^O fe 2 
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If additionally 



d 3 k 




(k)\ 2 


(27r) 3 2|fc| 




k 





< oo 



then Riemann-Lebesgue lemma (Chapter 2.5, Theorem 2.2 in 53]) implies 

d 3 k 

'O (k)\ 2 cos(kx)cos(\k\t ) = 
|0 (fe)| 2 sin(A;x)sin(|fc|to) = 0. 



Hence 



lim , , , , , 

t ^±oo J (27r) 3 |fc| 2 

d 3 k 

lim 



t ^±oo J (27r) 3 |fc| 2 



(565) 



lim (il>,N\</> c (t-t ,x,N)\ii>,N) = 

to— >±oa 



(566) 



for all states \ip,N) generated from vacuum by unitary transformations constructed from 
the iV-representation of HOLA. In this weak sense one can write 



lim <j) c (t-t ,x,N) = 0. 

to— >±oo 



(567) 



Let us note that asymptotic vanishing of compensating fields leads to the additional condition 
( 15651) that should be satisfied by vacuum wave functions. 

For comparison let us check what happens in the usual representation, Io(k) = 1, 

d 3 k 



(j) c (t-t ,x) = q 



■ cos(kx — \k\to) 



(2n) 3 \k\ 2 

Q 1 (2n) 3 \k\ 2 cos (l fc K*-M-fc-30 

poo p2tt ptt 

- 3 J dn J dip J sin # cos (n(t — t ) — \x\kcos8) 
Q 



o Jo Jo 

oo r—1 



(27T) 2 J A 



dn I du cos (^K,(t — t ) — \x\ku) 



(27T) 2 

g 

(2tt) 2 



I 

dn- — — sin (nit — t ) — lajl/tw) 
\x\k 



dn 



sin (k(£ — to + 


£c|)) + sin (/t(t - t + 






a? 


ft 



Note that I denote \k\ — K and dn, integration over the radial coordinate, which should not 
be confused with the measure dk on the light cone. For t = to 



</> c (0,a?) 



_q_ 

2vr 2 Jo 



dn 



sin k\ 


X 




K 


X 


Ati\x 
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For At = t — to 7^ the field is symmetric in time, c (At, x) = 4> c {— At, x). Assume At > 0, 



:(At,£c) 



(2vr) 2 



dn 



sin («(|cc| + At)) q 



\x\k 



+ 



(2vr) 2 Jo 



Q 


x\ 




At 


\x 


- At 


(2tt) 2 


\x 




\x 


-At 



(Ik 
dn 



sin [k(\x 


-At)) 


\x 




sin (k( \x\ 


-At)) 


k(\x\ - 


-At) 



? 1 



For \x\ = At > 



47r|cc| 2 



:(At,a;) 



(l + sgu(|x| - At)) = -r^-Mxl - At). 

\ / 47T \X\ 



(27T) 2 



sin2fi;kc q 7t 

dK : ; = r— = j r- 

k\x\ Aix 2 \x\ 2 An \x\ 2 



9 1 



For any t, to, 



c (t - t ,x) 



Ati\x 



■e(\x\ - \t-t \). 



For fixed t and \x\ ^ 



lim c (t — to, x) = 0. 

tg— >±oo 



As we can see, the step function is defined by 



for x < 



0(z) = < h for £ = • 



1 for x > 



(568) 



(569) 



(570) 



(571) 



(572) 



The role of compensating field is similar to the one of compensating currents occurring 
in some approaches to quantum electrodynamics (cf. 52j). A difference is that compen- 
sating currents are introduced in an ad hoc manner, whereas our compensating field is a 
consequence of quantum dynamics. 



B. Generalized Coulomb field 

In order to understand why I call <p g c{x) a generalized Coulomb field let us note that 
setting Io{k) = 1 one indeed arrives at the ordinary Coulomb field 

/d 3 k q 
(27r )3| fc |2 cos k ■ x = ~4^| = M*) = -0c (0, x). (573) 
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However, we know that Io(k) = 1 implies |Oo(fc)| 2 = 1 (since |Oo(fc)| 2 is the vacuum average 
of Io(k)) which would lead to infrared divergences and all types of internal inconsistencies 
of the formalism. So, take 

J (fc, N) = ^ (io(fc, 1) ® J (l) ® • • • ® /o(l) + ■ ■ ■ + /o(l) ® • • • ® 4(1) ® Io(fe, I! 

gC (a3, N) = — (4> gC (x, 1) <g> J (l) <8> • • • ® / (1) + • • • + Io(l) ® • • • ® io(l) ® 0gc(a3, 1) 

/d^ Jv f cos ife * 

(27T) 3 |fc| 2 C ° S(fc ' ^ ' fe ^ fe ' ® 1 = ~ 2<1 J ^ \k\ |fc ^ fc| ® 1 

„ cos • X fr .„A\ 

= -2q <g> 1, (574) 

|fe| 

k = J dkk\k)(k\. (575) 

Of course, <j) g c(x,N) is a field operator commuting with all elements of HOLA (the same 



modification of Coulomb's law was found in full electrodynamics in 27]). Vacuum averages 
of <j) g c{x, N) are ordinary functions that can be compared with 4>c{ x , N). 
Let us thus consider the average evaluated in a vacuum state \Oq, N), 

/d 3 k 
\O (k)\ 2 cos k-x. (576) 

Probability |O (fc)| 2 occurs in place where in the usual approach one might, if needed, put 
a cutoff function x(fc). Cutoff functions are assumed to vanish for large and small k but 
somewhere in between infrared and ultraviolet regimes achieve their maximal value 1. It is 
obvious, however, that the maximal value of |0o(fc)| 2 is greater from 0, but may be different 
from unity. It follows that |Oo(fc)| 2 itself cannot be regarded as the usual cutoff function. 
The latter is obtained if one divides |Oo(fc)| 2 by its maximal value Z, say. 

We know that Poincare transformations transform |Oo(fe)| into |Oo(£ -1 fc)|. Since the 
maximum 

Z = max{|O (fc)| 2 } = max{|Oo(-& _1 fc)| 2 } (577) 

k k 

is Poincare invariant we can define in Poincare covariant way the normalized "cutoff func- 
tion" 

X o(k) = \O (k)\ 2 /Z, 0< X o(fc)<l. (578) 

Accordingly, 

/d 3 k 
^ ^cosfc-a;, (579) 
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suggesting that qZ is the physical charge. 



studied in the literature so far (cf. 



20 



?his statement is almost true but all the examples 
271 ]) show that the correct physical parameter that 
should be compared with experiment is q 2 Z = (qy/Z) 2 = q 2 h . At this stage I propose the 
readers just to take this statement for granted, without proof. In other words, if we want to 
test the average with respect to the classical Coulomb field one should compare — q 2 h / (47r|aj|) 
with 

/d 3 k 
(27r)3|fc|2 x(fc) cosfc-x, (580) 

where g p h = q\fZ is the renormalized (i.e. physical) charge and Z is the renormalization 
constant (analogous to Z% known from textbooks [521]) . Let us also note here that condition 
( 15651) makes 

/d 3 k 
(27r)3 | fc | 2 |Qo(fc)| 2 (581) 

finite. 

To have a feel of general properties of q(O , N\(p gC (x, N)\O , N) let us take a rotationally 
invariant |Oo(fc)| 2 which vanishes for < \k\ < k\ and |fe| > k% but is otherwise constant 
(equal to some Z). Then Z can be then computed from 

l^/^W^^jT^^^M, (582) 
Z = (583) 



= gVZ = q-^=^. (584) 



Now 



q{O 0: N\(j) gC (x, N)\O , N) = -j^r dn d6 sin 6 cos(k cos 6\x\) 

(2tt) J kl J 



as; / ducos(/c aj u) 



4 rk2 



2tt 



sin k 


x\ 


K 


X 





9ph 


Si(A; 2 


cc|) -Si(fci 




Ati\x 


tt/2 



(585) 



The sine integral function (see Fig. 7) 
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0.5 




FIG. 7: Normalized sine integral ^Si(kx) for small x behaves as 2kx/7r but for large \x\ approaches 
1 independently of k (on the plot k = 10). This is why (Oq, N\c[)gp(x, N)\Oq, N) is finite at the 
origin and for large |x| decays faster than 1 / 1 ae | . If k\ were exactly (no infrared cutoff) then 
(Oq, N\(j) g c(x, N)\Oo, N) would asymptotically approach the Coulomb field. However, ki = 
would be in contradiction with lim^o Oo(k) = which we assume. 



Si(x) 



sin k 



dn 
o K 



has series expansion 



Si(z) = x-x 3 /18 + x 5 /600 -a; 7 /35280 + x 9 /3265920 + 



(586) 



(587) 



so 



lim Si(fc2X) " Si(fclX) = k 2 - kl . 
Another important property of Si(x) is the smallest solution of 

Si(oo) - Si(x) = 



(588) 



(589) 



which is numerically found to be x ~ 1.92645. What this physically means is that for a very 
large k 2 (a property one expects on physical grounds) the potential 

Att\x\ tt/2 1 ' 

changes sign if |cc| ~ 1.92645/A;i. Oscillating changes of sign are implied by oscillations 
of Si(x) depicted at Fig. 7. It is interesting that suggestions of analogous changes of sign 



have been serious 
Newton law (see 



y considered in cosmology and astrophysics in the context of gravitational 



54j for a systematized guide over the literature of the subject). 
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Current experimental data show that exact Coulomb law is indistinguishable, up to |cc| 
as large as the Earth-Sun distance, from the Yukawa law 

_ j£_ e -M/A (591) 
4n\x\ 

where A is of the order of the Earth-Moon distance (more than 300000 km, corresponding 
;o the Compton wavelength of a particle whose mass is 10~ 48 g < m < 10~ 47 g, cf. the review 
551] ). The experiments are supposed to test the value of photon's rest mass — the tacit 
assumption being that for massless photons and pointlike sources one necessarily gets the 
Coulomb law. However, we have just seen that the result depends also on quantization 
procedures — in my formalism the field is massless and the source is pointlike, and yet the 
Coulomb law is derived in a generalized form. 

Anyway, returning to the data, the Yukawa potential (I59ip effectively deviates from the 
Coulomb law for large \x\, a behavior controlled in my formula by the infrared cutoff k\. 
The data from [55] show that varying A between, roughly, A m ; n = 300000 km and infinity 
we remain within the experimental uncertainty even for |x| as large as some 10 9 km. The 
error of determining the exact Coulomb law thus can be defined as 

A(x) = r-r- — . (592) 
\x\ \x\ 

The scale of infrared cutoff is determined by k\, which can be estimated on the basis of 

<! ' HA " - AO.) = 2 eH " /Vi " - JL < ISiW-SW) = J_ _ 2 < J_ 

|x| \x\ \x\ \x\ 7r/2 |a?| 7r|cc| |cc| 

at least for \x\ < 10 9 km. So, 

g— |x|/A m j n 1 2 

2 — 2— < —SUUxl) < 0, 

\x\ \x\ n\x\ 

and finally 

< tt(1 - e-'*'/^) - Si(Jfei|a;|). (593) 

Fig. 9 shows that for k\ = 2/X min (or smaller) inequality H593I) is satisfied for all values of 
\x\. 

The conclusion is that even the most precise tests of the Coulomb law available so far 
do not contradict the possibility that Io(k) ^ 1. Oo(fe) may be a nontrivial function that 
tends to with k — > 0. The analysis I have presented is based on a concrete form of Oo(k) 
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£2=10,^=0.1 £2=10,^=0.01 




0123456 0123456 



£ 2 =100, ^=0.1 £2=100,^=0.01 




1 rrr-r-r-: -r-, 



0123456 0123456 

FIG. 8: Comparison of 1/x with 2[Si(/c2x) — Si(A;ix)]/(7rx) for arbitrarily chosen ki and hi- 

but one should not expect drastically different predictions if Oo(k) decayed in the infrared 
regime in a smoother way. 

Another important test of the Coulomb law would be the energy spectrum of hydrogen- 
like atoms with the operator potential (f) g c(x). I hope to return to this problem later on in 
these notes. 

C. "Photon statistics" 

Averages of generalized Coulomb field did not depend on the parameter N of the repre- 
sentation of HOLA. N becomes essential for the n- "scalar photon" statistics. 
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*i=6/A min 
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Let us begin with the unitary operator 

e -i H o dr V(r) = exp f_ iq t dr ^ Q iy) 
^ Jto 

= exp (-iqj^drj dk(c{k, N)e~ iWT + c{k, iV) t e i|fe|r 

-i, y ^(^jv) _. |fc| +C (fc,jv)t . |fc| 

= exp y dk(at )t0 {k)c{k, N) j - a ttto (k)c(k, N)j 

= V (a ttto ,N). (594) 
This is a coherent-state displacement operator with 

gi|fe|*o _ gi|fe|t 

= 9 j^] • ( 595 ) 

Repeating calculations analogous to those that led to (11841) we obtain statistics of "scalar 
photons" produced by a pointlike charge, with vacuum initial condition at t = to, 

p(n,N) = (O ,N\exp(- J dk\a tjto (k)\ 2 I (k, iV)) ^ ( J dk\ | a t)to (fex ) | 2 J (fei , A^)^ " | O , iV) 

= ^^(O ,iV|exp(A y dA;| aM0 (fc)| 2 J Q (fc,iV)J|O Q ,iV) ^ (596) 

Generating function occurring in (I596P can be further transformed, leading to Kolmogorov- 
Nagumo average 

(O ,iV|exp(A y ^|a M0 (fc)| 2 J (fc,iV))|O ,iV) = ( J dA;|0 (fc)|V^oWI 2 ) "(597) 
The limiting case 

lmi^ ( y dA;|Oo(fc)| 2 e A ^ |a «o( fe )l 2 )' V = exp (a y dk\O (k)\ 2 \a tAi (k)\ 2 ^ (598) 

is a regularized Poisson distribution. Explicitly, 

2 2 — e <[*[C*o— *) _ e -<[*[C*o— *) 2 2-2cos(|fc|(t-t )) 



l«Mo( fc )l = Q = ^ 



|fc| 2 * |fep 



2 sin 2 (|fc|(t-t )/2) 



= q 2 ^—L (599) 

(|fc|/2) 2 

and 

exp (A y dk\O (k)\ 2 \a t , t0 (k)\ 2 ) = exp^y ^fc | O (fc) | 2 ^ l° )/2) j 

A 2 [„ ^, sin 2 (|fc|(t-t )/2) 
= exp ^Ag ph y rffcx(fc) (|fc|/2) 2 ' 
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Let us take a closer look at the exponent in (I600p . Since, 

sin 2 (|fc|At/2) A 9 sin 2 (|fc|At/2) A „ . , 

< 77- 1 , J = At 2 . ' ' < At 2 601 

" (|fc|/2)2 (\k\At/2y ~ y } 

the integral is finite at least for finite At = t — to, 

s ^ |Q ° (fc)| >2( y /2> - {t ' to)2 i dm[k)? = {t ~ to)2 - (6 ° 2) 

In the ordinary approach we would obtain the same formula (16001) but with |Oo(fc)| 2 = 1, 



i.e. 



d 3 k sin 2 (|fc|At/2) 16tt f°° , sin 2 (KAt/2) ,„ rtn , 

aK . (603) 



(2tt) 3 2|A;| (|fc|/2) 2 2(2vr) 3 J k 

The latter integral vanishes for At = 0, but becomes divergent for At > 0. More precisely, 
for < ki < fc 2 , 

fe2 J sin 2 K Ci(2A; 1 )-lnA; 1 -Ci(2A; 2 ) + In k 2 , cnA . 

clk = . (604) 

fci K 2 

The infrared limit is finite, 

lim (Ci(2fci) - lnfci) = 7 + In 2 (605) 

ki->0 v ' 

(Ci and 7 ~ 0.577216 are, respectively, the cosine integral function and the Euler constant). 
The infinity comes from the ultraviolet divergency of ln&2- An infrared divergency would, 
however, also appear if the charge was accelerating (for details of standard calculations 



cf. [52[ and Sec. 4-1-2 in 56]; analogous results for iV-representations of HOLA and true 
electromagnetic fields are described in detail in |27|). One concludes that the standard 
approach to pointlike classical sources leads to mathematically ill defined objects such as 
e Ao ° and thus is mathematically inconsistent. 

In "my" formalism one finds that probability of finding n "scalar photons" is given by 
the vacuum average of POVM 

U(n,N) = exp(-y dk\a t>to (k)\ 2 I (k, jV)) ^ ( J dfciKJfci)! 2 ! / (fci, N)J 

1 d n / f \ 
= ^^T ex P ( A / dk\a t!to (k)\ 2 I (k,N)j a= _^ (606) 

00 

J>(n,iV) = J (JV). (607) 

n=0 
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This POVM is uniquely and well denned, but its average depends on the choice of vacuum, 
of course. 

Now let us discuss the issue of the limit t — > ±00. Let us return to 

,2 -2cos (\k\{t-t )) 



dk\0 (k)\ 2 \a t , to (k)\ 2 = q 2 j dk\O (k)f 



\k\ 



2q 2 



d 3 k 


\Oo(k)\ 2 


(2vr) 3 2|fc| 




k 


2 



_2g 2 / d3ft |Q (jfc)l2 COS (l fc l*) COS (l fc l*0 



(27r) 3 2|fc| 1 uv " |fc| 2 
I (^|Oo(^)l 2 5m(|fc|i |i S | i " ( "" f ° > « 



If the first integral 



d 3 k 


\Oo(k)\ 2 


(27r) 3 2|fc| 




k 


2 



is finite (compare (14421) ) then |Oo(fc)| 2 /|fc| 3 satisfies assumptions of the Riemann-Lebesgue 
lemma, and 



lim / dk\O (k)\ 2 \a t , (k)\ 2 = 2q 2 [ 



d 3 k 


\O (k)\ 2 


(27r) 3 2|fc| 




k 
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